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^si ' Abstract. The eigenvalue statistics of a pair (Mi, Af2) of n x n Her- 

mitian matrices taken random with respect to the measure 

-^ exp ( - nTr(V'(Afi) + WiAh) - rMiMa)) dA/jdMa 

("V^ ' can be described in terms of two famihes of biorthogonal polynomials. 

C^ I In this paper we give a steepest descent analysis of a 4 x 4 matrix- valued 

Riemann-Hilbert problem characterizing one of the families of biorthog- 
' (-h '' onal polynomials in the special case W{y) — j/*/4 and V an even poly- 

Oj, nomial. As a result we obtain the limiting behavior of the correlation 

kernel associated to the eigenvalues of Mi (when averaged over M2) in 
the global and local regime as n — » 00 in the one-cut regular case. A 
C^ , special feature in the analysis is the introduction of a vector equilibrium 

problem involving both an external field and an upper constraint. 
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1. Introduction 

1.1. Two matrix model. The two matrix model in random matrix theory 
is a probability measure 

(1.1) — e-"^ {v{M,)+wiM2)-rM,Ah) dMidMa 

defined on pairs (Mi, M2) of n x n Hermitian matrices. Here Z^ is a normal- 
ization constant, V and W are two polynomials of even degree and positive 
leading coefficients, r 7^ is a coupling constant and dMj denotes the prod- 
uct of the Lebesgue measures on the independent entries of Mi for i = 1, 2. 
The model was introduced in [49^ [64] as a generalization of the unitary one 
matrix model which allows for a larger class of critical phenomena, see also 

2 
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[5T] and [371 EH] for a survey. For more recent advances in the physics 
Uterature, see e.g. [U HSj [45], [46] and the references cited therein. 

An open problem in random matrix theory is to give rigorous asymptotic 
results on the eigenvalue statistics of Mi and M2 as n — > (X). A natural 
approach is to use the connection with biorthogonal polynomials. In [66] 
it is shown that the eigenvalue statistics can be described in terms of two 
families of polynomials {pk,n)k and {qi^n)i, where pk^n and qi^n are monic 
polynomials of degrees k and /, respectively, satisfying 



(1.2) // PM(^)'?^,n(y)e-"(''(^')+^(^)-™'^) 



dx dy = 0, k ^ I. 

These polynomials are well-defined and have real and simple zeros [5T] . From 
a complete asymptotic description of the polynomials, it is possible to com- 
pute the limiting behavior of the eigenvalue statistics. This has been carried 
out for the case 

(1.3) V{x) = x'^/2 and W{y) = ay^/2 

in [H] . Although heuristic calculations can be found in the physics literature 
[42j . fully rigorous asymptotic results for the biorthogonal polynomials for 
more general V and W are not known. 

The orthogonal polynomial approach to random matrices proved to be 
successful in the one matrix models. The orthogonal polynomials appearing 
in these models are characterized by a 2 x 2 matrix valued Riemann-Hilbert 
problem [47]. In [33l [34] the authors applied the Deift/Zhou steepest de- 
scent method to this Riemann-Hilbert problem and obtained a complete 
asymptotic description of the polynomials. As a result it was possible to 
give a rigorous proof of the universality conjecture for the local eigenvalue 
correlations. 

Inspired by the success of [33l [34] , several attempts were made to study the 
asymptotic behavior of the biorthogonal polynomials by Riemann-Hilbert 
methods. A number of Riemann-Hilbert problems for the biorthogonal poly- 
nomials were formulated [HI [HI [50l [56] , but a successful steepest descent 
analysis has not been carried out so far. It is included in Deift's list of 
major open problems in random matrix theory and the theory of integrable 
systems [31] . 

In this paper we present the first complete steepest descent analysis for a 
case beyond the fully quadratic case (jl.3p . We analyze the Riemann-Hilbert 
problem for the biorthogonal polynomials pn^n given in [56| for the special 
case 

(1.4) W{y) = y /4 and V is an even polynomial. 

As a result we are able to compute asymptotics of the eigenvalue correlations 
of the matrix Mi , when averaged over M2 ■ 
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1.2. Unitary ensembles. Let us first recall some aspects of the unitary 
ensembles and orthogonal polynomials. In the unitary ensemble one consid- 
ers nxn Hermitian matrices taken randomly with respect to the probability 
measure defined by 

(1.5) ^e-"^^WdM, 



where V is such that 



X 



lim - — -^ — = -l-cx). 

-^±oo log(j;^ + Ij 



Let pk^n be the unique monic polynomial of degree k satisfying 

/•oo 

(1.6) / pfc,„(x)x^e-"^(") dx = 0, j = 0,...,k-l. 



— oo 



Then the eigenvalues of M describe a determinantal point process on M with 
kernel Kn defined by 

n 7^ T^ f \- 2 _„^W+n«) Pn,n{^)Pn~l,n{y) " Pn,n{y)Pn-l,n{^) 
[I. I) J\n[X,y) — Jn-1^ ' 

x-y 

where the constant jn-i is the leading coefficient of the orthonormal polyno- 
mial of degree n — 1. Thus the joint probability density V on the eigenvalues 
xi, . . . ,Xn is (up to a constant) equal to the determinant 

(1.8) V{xi,... ,Xn) = — det(K„(xj,Xj))"^.^-^, 

and similarly for the marginal densities for k = 1, . . . ,n — 1, 



(1.9) 



(n - A;)! u 

V{xi,...,Xn) dxk+i---dxn = j — ^ det (i^„(xi, x,' )).'._-, . 

n—k times 

In order to compute the asymptotic behavior of the eigenvalue statistics, it 
is sufficient to obtain asymptotics for the orthogonal polynomials and the 
kernel Kn- 

Fokas, Its, and Kitaev [47] characterized the orthogonal polynomials p„^„ 
in terms of a Riemann-Hilbert problem (RH problem) . It consists of seeking 
a 2 X 2 matrix valued function Y satisfying 

Y is analytic in C \ M, 

y+(x)=y_(x)(^J " ^ j, xGM, 

y(z) = (/ + 0(z"i))("J ^°„), z^oo. 

Here Y+ and Y- denote the limiting values of y on R when M is approached 
from above and below, respectively. The unique solution of the RH problem 



(1.10) 
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is given by 

/-, -|-|N Y= ( ^"-'^ C [pn,n^ " j 

where C f denotes the Cauchy transform 

(1.12) [Cf){z) = — r -I^dx, zeC\R. 

In [23] and [33], the Deift/Zhou steepest descent method for RH problems 
[36] was apphed to obtain the asymptotic behavior of 1" as n ^ oo, and 
hence of the orthogonal polynomials. 

A key ingredient in the steepest descent analysis in [33l [3l] is an equilib- 
rium measure, see also [301171]. This measure is the unique minimizer of the 
energy functional 

(1.13) If log -^ d/i(x)d^(y) + [ V{x)dfi{x) 
JJ \x-y\ J 

among all Borel probability measures fi on M. If V is real analytic then 
the equilibrium measure is supported on a finite number of intervals, has an 
analytic density in the interior of each interval and vanishes at the endpoints 
[32] . It is the weak limit of the normalized counting measure on the zeros 
of Pn,n as n — >■ oo. Moreover, it describes the limiting mean eigenvalue 
distribution of a matrix from (jl.5p . 

The papers [331 IM] had a major impact on the theory of random matri- 
ces and orthogonal polynomials. Inspired by these papers, several authors 
extended the methods of [331 El] to obtain asymptotics for different types 
of orthogonal polynomials. For example, for orthogonal polynomials on the 
half line [74] , on the interval [57] , and on the unit circle [621 [63] . Another 
important development is the asymptotic analysis for discrete orthogonal 
polynomials ^. In all these cases the orthogonal polynomials can be char- 
acterized in terms of a 2 x 2 matrix-valued RH problem and an associated 
equilibrium measure plays an important role. 

1.3. Two matrix models and biorthogonal polynomials. Let us now 

return to the two matrix model (jl.ip and the biorthogonal polynomials pk^n 
and qi^n in (|1.2p . In the one-matrix model the eigenvalues of the random 
matrix follow a determinantal point process on M whose kernel is the re- 
producing kernel corresponding to the orthogonal polynomials. In the two 
matrix model we have a similar result but the situation is more complicated. 
Define the transformed functions 

(1.14) Qk,n{^) = e-"^(-) r (?M(y)e""('^^'^"™')dy, 

(1.15) Pk,n{y) = e-"^(^^ r PM(^)e""(''^"^"""')dx, 
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and let h^. ^ be defined as 

(1.16) hl^ = r r PM(^)9fc,n(2/)e-"(^(^)+^(^)-™^) dxdy. 

J—oo J —oo 

Consider the following kernels 

n— 1 ^ 

(1.17) K|"^(xi,X2) = ^ T2— Pfc,n(a;i)Qfc,n(a;2), 

k=0 ^'"i^ 
n—1 ^ 

(1.18) K^2 (yi^y^) = "^ j^PkAyihkAvi), 

n—1 ^ 

(1.19) K[2\x,y) = '^j^Pk,ni^)lk,n{y)^ 

k=0 '''"' 

(1.20) Ki-\y, x) = J2 T^PkAy)QU^) - e-"(^(-)+^(^)— ^) . 

fc=0 '^.'^ 
Then the joint probability density function V for the eigenvalues xi, . . . ,Xr^ 
of Ml and the eigenvalues yi , . . . , y^ of M2 is given by [121 [66] (see also [65 
Chapter 23]), 



(1.21) 'P{xi,...,Xn,yi,...,yn) 



1 



i>fi"^(xi,Xj))".^i (^S^(a;*,yj))r,=i' 



^^^ . .„., 



(n!)2 {{Kif{y,,x,))l^^^ {K^^ {y^y.^J ' 
Moreover, the marginal densities also have the determinantal form 

(1.22) /■■ ■/ ^(a:i,...,Xn,yi,...,yn)dxA;+i---dxndyi+i---dy„ 

n—k+n—l times 

^ (n-/)!(n-fc)! /(<)(x.,x,))J (^S^(-.,%)) ■ ' =A 

(-0^ \{<\y^,x,))^C-^ {4-\y.,y,))[^J ■ 

After averaging over the eigenvalues of M2 , we see that the eigenvalues of Mi 

(n) 

follow a determinantal point process with kernel Kj^^ . Similarly, the eigen- 

(n) 

values of M2 follow a determinantal point process with kernel -f^22 • Both 
these determinantal point processes are examples of biorthogonal ensembles 
in the sense of [21]. 

To determine the asymptotic behavior of the eigenvalues of Mi and M2 
it is sufficient to determine asymptotic behavior of the biorthogonal poly- 
nomials and the kernels -fC)" . 

The biorthogonal polynomials have been studied for many years and many 
interesting properties have been discovered [H [9l [lOl [TTl [12l [16l \^ HSJ 
[50l [56] . Although heuristic results on the asymptotic behavior of the 
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biorthogonal polynomials can be found in [32j, rigorous asymptotic results 
have not yet been obtained. 

A first step to an asymptotic analysis by means of RH methods is the 
formulation of a RH problem for the biorthogonal polynomials |121 |4H \50\ 
I56j . Here we follow |56j . The relationship with the RH problem of [12] has 
been clarified in [16]. In fact, this relationship will be exploited as part of 
our analysis. 

Assume W is a polynomial of degree d\Y and define dw — 1 weights 

(1.23) 



Wj^nix) 



-nV{x) 



yiQ-n{W{y)-Txy) ^^^ 



j = 0,l,...,d' 



w 



2. 



Then the RH problem associated with the biorthogonal polynomial p^^n is 
the following. We look for a dw x dpv'-matrix valued function Y satisfying 
the following properties 
(1.24) 

' Y is analytic in C \ M, 

1 ••• 



y+(x) = y_(x) 



Vo 







X G 



Y{z) = {l + 0{l/z))diag{z'',z-''\...,z 



oo. 



where kj is the integer part of {k + d\Y — 1 — j)/idw — !)• In [SS] it is proved 
that this RH problem has a unique solution given by 



/ 



(1.25) 



Y 



Pk,n 

(0) 

Pk-l,n 



> Jdw-2) 
\i^k-l,n 



C{pk,nWo^n) 
C{Pkll,n'^0,n) 



Cipit~"^^o, 



C{pk,nWd^-2,n) \ 
C{ptiwdy^-2,n) 



C{p^k-l,n^Wdw-2,n)) 



Aj) 



where pk^n is the monic biorthogonal polynomial of degree k and p^_i^, 
j = 0, . . . , dw — 2, are certain polynomials of degree < A: — 1. Here Cf 
denotes the Cauchy transform as given in (I1.12p . 

In the case dw = 2 the RH problem (ll.24p reduces to the 2x2 ma- 
trix valued RH problem for orthogonal polynomials. Then the biorthogonal 
polynomials pk,n are simply orthogonal polynomials on the real line with 
respect to a varying exponential weight, see also [41] . 

For dw > 2 the biorthogonal polynomials pk^n do not reduce to orthogonal 
polynomials, but instead they are examples of what is known as multiple 
orthogonal polynomials. Multiple orthogonal polynomials for r weights are 
characterized by (r + 1) x (r + 1) matrix valued RH problems [73], and (|1.24p 
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is an example of such a RH problem. The steepest descent analysis has not 
been applied to the RH problem ()1.24p . 

Certain other systems of multiple orthogonal polynomials and their asso- 
ciated RH problems were successfully analyzed recently. These polynomials 
were related to matrix models with external source [31 [T9l [20l [HI] and non- 
intersecting paths [27l [28l [54] . See also [HI [14] for a 3 x 3 matrix valued 
RH problem describing a Cauchy two matrix model. The steepest descent 
analysis of the larger size RH problems revealed several new features that 
are not present in the analysis of the 2x2 RH problem (jl.lOp . 

In our opinion, an important obstacle in the asymptotic analysis of ()1.24p 
is the lack of a tractable equilibrium problem. One of the main contribution 
of the present paper is a suitable equilibrium problem that we use in the 
steepest descent analysis of ()1.24p for the special case ()1.4p . 

1.4. RH problem for case W{y) = y^/4. In this paper, we will not treat 
the general two matrix model with polynomial potentials, but we restrict 
ourselves to the special case (jl.4p . Due to the assumption W{y) = y^/4 
the RH problem ()1.24p is of size 4x4. The assumption that V is an even 
polynomial introduces a symmetry with respect to the imaginary axis into 
the problem. This will be important for the steepest descent analysis. We 
also make the additional assumption that re is a multiple of three. This is 
not essential and is made only for reasons of exposition, since it simplifies 
many of the formulas. 

Due to these assumptions the RH problem ()1.24p characterizing the biorthog 
onal polynomial pn,n takes the form 

Y is analytic in C \ M, 

1 



(1.26) 



Y+{x) = Y^{x) 



/I 







Y{z) = il + 0{l/z)) 









Vo 





1 







-n/3 









^— n/3 





W2,nix)\ 




1 





-n/3l 



X G 



/ 



oo. 



with 

(1.27) 



w 



3,n 



(x) 



-nV{x) 



Jp-n{v*/4:-Txy) 



ye 



dy, 



0,1,2, 



This RH problem has a unique solution which is given by ()1.25p . A central 
observation in [56j that leads to this result, is that the polynomial pn,n is 
a multiple orthogonal polynomial with respect to the three weights Wj^n on 
the real line. This means that pn,n is the unique monic polynomial of degree 
re such that 



(1.28) 



Pn,n{x)x^Wj^n{x) dx = Q, 1 = 0,... 



re 
3 



1, j = 0,l,2. 
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The polynomials pk,n are multiple orthogonal polynomials of type II. 
There are also multiple orthogonal polynomials of type I. These appear 
in the function Qk,n of (|1.14p . which can be written as 

3 

(1-29) Qk,n{^) = 5]^g„)(x)t/;,-4x), 

3=0 

for certain polynomials A,V ^, called multiple orthogonal polynomials of 
type I. In |26] a Christoffel-Darboux formula is proved for the reproducing 
kernel associated to multiple orthogonal polynomials in a general setting. In 
our case it applies to the kernel i^j" of ()1.17p . As a result of the Christoffel- 

Darboux formula the kernel K\^ can be expressed in terms of the solution 
Y to the RH problem ()1.26p as follows. 

Proposition 1.1. We have that 
(1.30) 

/1\ 



,{n), . 1 



Kii{x,y) = 2^Ux_ ) (Q ^o,n(y) wi,niy) w2,niy))y+^iy)Y+{x) 

for x, y G M. 

Proof. See t26j. D 

2. Statement of results 

Our results deal with the limiting behavior as n ^ cxd of the eigenvalues 
of the matrix Mi in the two matrix model 



J_ -nTr {v{Mi)+W{M2)-tMiM2) 



Zn 



g-n ir yv (iviij + vv {ivi2)-Tivnivi2J Jj^^,^^^^ 



with T > under the following assumptions 

(a) y is an even polynomial, 

(b) Wiy) = 2/V4, and 

(c) n is a multiple of three. 

Guionnet [l8] showed that the limiting mean eigenvalue density of the 
matrix Mi exists. Her result is valid in much greater generality than stated 
in the next proposition. 

Proposition 2.1. The limiting mean density of the eigenvalues of the ma- 
trix Ml exists. That is, there exists a probability measure pi on M with 
respect to Lebesgue measure such that in the sense of weak convergence of 
measures 

(2.1) lim -k\'^\x,x) dx = d^i. 

n— >cx3 n 

Proof See [48j. D 
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Guionnet characterized ni by a variational problem. Under the above as- 
sumptions (a)-(b), we are going to characterize /xi in terms of an equilibrium 
problem from potential theory. It is not clear to us if the two minimization 
problems are related to each other. 

2.1. The equilibrium problem. As already mentioned, we characterize 
^1 by means of an equilibrium problem from potential theory. Main refer- 
ences for potential theory in the complex plane are [70^ [7T] . 

We will work with (non-negative) measures on C, but only on R and iM. 
A signed measure is real and can be negative. The support of a measure v is 
denoted by S^v). The measure may have unbounded support. If a measure 
1/ has unbounded support then we assume that 

(2.2) / log(l-h|2;|) di/(x) <oo. 
The logarithmic potential of u is the function 

(2.3) U^iz) = /log I ^ A i^ix). 
The logarithmic energy of v is defined as 

(2.4) i^^)= ffiog^-^diy{x)du{y) 

J J \x-y\ 

and I{i') G (— oo, -|-cxd]. If z^i and z^2 are two measures with finite logarithmic 
energy, then their mutual logarithmic energy is defined by 

(2.5) I{i^i,i^2)= tt\og—^^dyi{x)dv2{y). 

J J \x-y\ 

For a given V and r > 0, we study the energy functional Ey defined by 

(2.6) Ev{yi,V2,v^) 

3 2 



^I(z.,)-^/(i/„z.,+i)+ j (v{x) 

.7 = 1 .7 = 1 '' ^ 



Y'M^'^ ) d^^i(^) 
i=i r- ' 



where i^i, v^^ i^3 are three measures with finite logarithmic energy. The 
equilibrium problem is the following. 

Definition 2.2. The equilihrium problem is to minimize Eylui, z/2, fs) among 
all measures vi, vi, vj, such that 

(a) the measures have finite logarithmic energy; 

(b) z/i is a measure on M with i^i(M) = 1; 

(c) z/2 is a measure on iM with 1^20^) = 2/3; 

(d) z/3 is a measure on M with f^iM) = 1/3; 

(e) z/2 satisfies the upper constraint 1^2 <^ o" where a is the (unbounded) 
measure on iM defined by 

(2.7) daiz) = ^T^/^\z\^/'^\dz\, zGiM, 
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where \dz\ denotes the arclength on iM. 

An electrostatic interpretation of the equilibrium problem for the energy 
functional ()2.6p is the following. Consider three types of charged particles. 
The particles of the first type are put on M and have total charge 1. The 
particles of the second type are put on iM and have total charge 2/3. The 
particles of the third type are put on R and have total charge 1/3. Parti- 
cles of the same type repel each other. The particles of the first and the 
second types attract each other with a strength that is half the strength 
of the repulsion of particles of the same type. So do the particles of the 
second and the third types. Particles of the first and the third types do not 
interact directly. Particles of the first type are influenced by an external 
field depending on V and r. 

In the equilibrium problem the particles distribute themselves in order to 
minimize their energy under the extra condition that the particle density 
of the second type particles does not exceed the density of a, where a is 
the given measure (j2.7p . Thus a acts as an upper constraint on the second 
measure. 

Equilibrium problems for a vector of measures with mutual interaction as 
in (j2.6p arise for Nikishin systems [68] in the theory of rational approxima- 
tion, see also the survey [2j, and the more recent papers [H [6l HO t [54 t [60 l [73] . 

Equilibrium problems with constraint appeared before in asymptotic re- 
sults for discrete orthogonal polynomials (Sj El EU [SU [69] , singular limits 
of integrable systems [351 , and convergence results for Krylov methods in 
numerical linear algebra j52j . 

We prove the following. 

Theorem 2.3. Let V be an even polynomial and r > 0. Then there is a 
unique niinimizer (fii, fj,2, fJ-s) of Ey subject to the conditions (a)-(e) in the 
equilibrium problem. Moreover, 

(a) The measure fii is supported on a finite number of disjoint intervals 
Uj=i[oi; ^j] with a density of the form 

(2.8) ^^ = hj {x)^{bj -x){x-aj), x G [oj , bj], 

where hj is real analytic and non-negative on [oj, bj], for j = 1, . . . ,N. 

(b) S{p2) = ilR CLiT'd there exists a constant c > such that 

S{a — ^2) = iM \ (— ic, ic). 

The measure o" — ^2 has an analytic density on iM \ {—ic, ic) which 
vanishes as a square root at ibic. Moreover, the logarithmic potential 
[7^2 {g such that 
q 

(2.9) -r'^/^|x|''/^ + U^^{x) is real analytic on M. 

(c) S{iJ,3) = M and 113 has a density which is analytic in M \ {0}. 
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(d) all three measures jii, ^2 oind fj,^ are symmetric in the sense that 
fij{—A) = Hj{A) for j = 1, 2, 3, and for every Borel set A. 

The symmetry property in part (d) of Theorem l2.3l is a direct consequence 
of the uniqueness of the minimizer and the fact that V is even. We state 
it exphcitly here, since it wih be used many times in what foUows, often 
without exphcit mentioning it. 

For given V2, v^, we can minimize EY{vi,V2^i^'i) with respect to vi only. 
Then we look for the minimizer of the energy functional 

(2.10) I{vi) + / {v{x) - ^r^/3|a;|4/3 - U^-^{x)\ dui{x), 

among all probability measures vi on M. This is a usual equilibrium problem 
with external field [301 EI]- The external field is analytic in M \ {0}, but 
possibly not at 0. However, due to (|2.9p the external field is analytic at if 
J^2 = /U2 is part of the minimizer (/^i, /i2, /U3) of the full equilibrium problem. 
Thus ^i is the minimizer of an energy functional with analytic external field 
on M. The statements in Theorem 12.31 about the structure of the measure 
/ii then follow from results of [32]. 

If we minimize only with respect to 1/3, with ui and V2 fixed, then we 
minimize the energy functional 

(2.11) I{v^)- ju^^{x)dv^{x), 

among all z/3 on M with total mass 1/3. Again this is an equilibrium problem 
with external field, but the external field —U^'^{x) is only slowly growing at 
infinity. In [72] such an external field is called weakly admissible. It is a 
consequence of the slow growth that the support 5(1/3) of the minimizer will 
be unbounded, and in fact it will be the full real line. 

Finally, if we fix ui and 1/3, and minimize with respect to ^2, then we look 
for the minimizer of the energy functional 

(2.12) I{v2) - I iU'^'ix) + C/'^3(x)) di/2(x), 

among all z/2 on iM with total mass 2/3 satisfying the constraint 1/2 < c. Here 
we again have a weakly admissible external field, but in addition there is 
the upper constraint. Again, we will have that the support of the minimizer 
is unbounded, and indeed it is the full imaginary axis. The constraint a is 
active on a symmetric interval [— ic, ic] around the origin. 

2.2. Variational conditions and regular/singular behavior. From the 
discussion above we see that each of //i, /X2, /U3 is the minimizer for an 
equilibrium problem with external field and/or upper constraint, and as 
such they are characterized by the following set of variational conditions. 

Proposition 2.4. The measures /ii, /X2 and /X3 satisfy 

(2.13) 2U^''{x) = Uf'^x)-V{x) + ^T^/^\x\^/^+£, x £ S{fii), 
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(2.14) 


2^^1(2;) > U^'^x) - Vix) + |t^/3|x|4/3 +^, 


for some 


constant i, 


(2.15) 


2C/^2 (x) = W^ (x) + W' (x), X e 5(a - ^2 


(2.16) 


2?7^2 (^) < ^^1 (2;) + f//.3 (^)^ X G iM \ 5(f7 ■ 


and 




(2.17) 


2[/^3 (X) = [/'^^ (^)^ ^ g 5(^3) = ^, 



x£R\S{fii), 



^2), 



Note that the inequahty in (j2.16p is strict. By contrast the inequahty in 
()2.14p need not be strict. Then as in the case of the one-matrix model and 
orthogonal polynomials |33j we can distinguish between regular and singular 
cases as follows. 

Definition 2.5. The measure /Ui is called regular (otherwise singular) if the 
functions hj in (|2.8p satisfy hj > on [aj, bj] for every j = 1, . . . , N, and if 
the inequality in (j2.14p is strict. 

A point X* E M ai which one of these conditions fails is a singular point. 
There are three types of singular points, namely 

• X* is a singular end point if x* £ {oj, bj} for some j = 1, . . . ,N and 
hj{x*) =0, 

• X* is a singular interior point if x* G {oj, bj) for some j = 1, . . . ,N 
and hj{x*) = 0, 

• X* is a singular exterior point if x* G M \ S{fii) and equality holds 
for X = X* in ()2.14p . 

We see that the structure for the measure /Ui in the global regime is the 
same as what is known for the one-matrix model in a real analytic external 
field [33]. It is known that the regular case holds generically in one-matrix 
models [55] and all three types of singular behavior can occur. The same 
is true in the present two matrix model, but the maybe surprising fact is 
that no other type of singular behavior can occur. We emphasize that this 
is related to our assumption that W{y) = y^/4. 

We conclude with an easy to check convexity result, which is analogous 
to the well-known fact that for a convex external field the equilibrium mea- 
sure is supported on one interval. The following is a similar result for the 
equilibrium problem of Definition 12.21 

Proposition 2.6. Suppose that V is even and that x 1— > V{y/x) is convex 
for X > 0. Then S{fii) consists of one or two intervals. 

The proposition applies to all even quartic potentials V{x) = ax^ -|- 6x^ 
with a > 0. 

2.3. One-cut regular case. The core of the present paper is a Deift/Zhou 
steepest descent analysis for the RH problem for the biorthogonal polyno- 
mials. As a particular result we are able to obtain the limiting behavior of 
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the kernel Kj^ as n — > oo. The steepest descent analysis simplifies in the 
one-cut regular case, that is, the measure ni is regular and supported on 
one interval. We will only deal with this case in this paper, but our methods 
can be extended to the non-regular and multi-cut situations. 

The non-regular cases for unitary ensembles have been treated recently, 
see |18t [23| ^^ for singular interior points, [25] for singular endpoints and 
[151 [22| [67] for singular exterior points. All three singular cases can appear 
in our situation, and can be implemented within our methods to obtain 

(n) 

the limiting behavior of the polynomials pn.n and the kernel K^^ for the 
singular cases as well. We will not go into details here. 

As part of our analysis we will express the minimizer (/.ii, /i2, ^3) in terms 
of a meromorphic functions defined on a four sheeted Riemann surface. 
The genus of the Riemann surface equals A^ — 1 where N is the number 
of intervals in the support of fii. Hence the genus is non-zero in the multi- 
cut case, which complicates the construction of the outside parametrix in 
the steepest descent analysis. For reasons of exposition we only deal here 
with the one-cut case and plan to return to the multi-cut case in future 
work. 



2.4. Global eigenvalue regime. Our first main theorem states that the 
measure fii that we obtain as part of the minimizer {ni, fj,2, fJ'3) for the equi- 
librium problem for Ey is the limiting mean distribution for the eigenvalues 
of Ml. We prove this in this paper for the one-cut regular case. 

Theorem 2.7. Let W{y) = y^/4, let V be an even polynomial and r > 0. 
Let (/ii,/i2)/^3) be the minimizer of the equilibrium problem described in 
Subsection \2.1\ above. Suppose that the measure /ii is one-cut regular. 
Then we have as n — > c« with n = 0(mod3) 

(2.18) lim -K-f")(x,x) = ^(x), 

n-^co n ax 

uniformly for x G M. 

2.5. Local eigenvalue regime. In the local eigenvalue regime we obtain 
the same universal limiting behaviors that are known from one-matrix mod- 
els. This is in agreement with the universality conjecture in random matrix 
theory, which says that local eigenvalue correlations in random matrix mod- 
els with unitary symmetry do not depend on the particular features of the 
model, but only on the global regime. 

We again restrict ourselves to the one-cut regular case. So we assume that 
the measure /ii from the solution of the equilibrium problem is supported 
on one interval and we write S{fj,i) = [—a, a]. 

Theorem 2.8. Let W{y) = y^/4, let V be an even polynomial and r > 0. 
Assume that the measure fii is one-cut regular. 
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(a) Let X* E (—a, a) and define 

(2.19) p:=^(^*)>o. 

ax 

Then, as n —> cxd with n = 0(mod3) we have for every /c G N that 



(2.20) lim det ( -^k['1^ ( x* + ^, x* + ^ 



k 



n^oo \ pn \ pn pnjj^-^^ 

sin7r(nj — Uj)) 



det 



Tl{Ui 



i,j=i 



uniformly for (ui, . . . , Uk) in compact subsets ofM.^. 
(b) Let p > be such that 

^^^ = ^(a-x)i/2(i + 0(a-x)), asx/a. 

ax vr 

Then, as n —> oo wii/i n = 0(mod3) we have for every k £N that 
(2-21) Jim det (i^< (« + ^,„+ ^))' 



det 



Ai{ui) Ai' (uj) - Ai'('Uj) Ai('Uj; 

uniformly for (ui, . . . , u^) in compact subsets ofM.^ . Here Ai denotes 
the usual Airy function. 

In Theorem 12.81 we only deal with the one-cut regular case. The non- 
regular cases that can occur in our situation, also appear in the unitary 
ensembles. In these cases, we also obtain the same limits for the kernel Kj^' 
as in the unitary ensembles. For example, if the density of the measure 
pi vanishes quadratically at the an interior point of the support, then the 
limiting behavior of the kernel K)^, in a double scaling limit, is expressed 
in terms of the ^-functions associated to the Hastings-McLeod solution of 
the Painleve II equation, see also [T8|[23]. 



2.6. Overvie'w of the rest of the paper. The proof of Theorem 12.81 fol- 
lows from a Deift/Zhou steepest descent analysis applied to the RH problem 
for biorthogonal polynomials. The Deift/Zhou steeepest descent analysis 
consists of a sequence of explicit and invertible transformations 

(2.22) Y^X^U^T^S^R. 

The RH problem for Y is not in a suitable form for an immediate appli- 
cation of the equilibrium measures from the equilibrium problem. In the 
first transformation Y >-^ X we transform the RH problem for y to a RH 
problem for X that is more suitable for further asymptotic analysis. This 
transformation depends on a method that we learned from the authors of 
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[16] . The construction involves Pearcey integrals and is presented in Section 
3. 

In the second transformation X i— > [/ we will use the equilibrium measures 
/ii, fi2 and /U3 from the equilibrium problem of Definition 12.21 and their 
corresponding y-functions. In Section 4 we discuss the equilibrium problem 
and we prove Theorem 12.31 and Propositions 12.41 and 12.61 

The properties of the gf-functions are conveniently expressed in terms of 
functions that come from a four sheeted Riemann surface, which in the one- 
cut case has genus zero. We introduce the Riemann surface in Section 5. 
The second transformation of the RH problem X ^^ U is given in Section 6. 

In the transformation [7 1— > T 1— > 5 in Section 7 we open lenses. In the 
transformation f7 1— > T we open the lens around the supports of 1^2 and /is. 
These supports are unbounded so special care has to be taken at infinity. In 
the transformation T 1— > S we open the lenses around the support of fii. 

The next step is the construction of a parametrices in Section 8. Here 
the Riemann surface is used again. We give a rational parametrization 
which allows us to give explicit formulas for the outside parametrix. In the 
multi-cut case the Riemann surface has higher genus and the construction 
of the outside parametrix is more complicated. We construct local Airy 
parametrices around the branch points of the Riemann surface. These are 
the endpoints ita of the support of /xi and the endpoints itic of the support 
of o" — /i2- The final transformation S >-^ R is also given in Section 8. It 
leads to a RH problem for R with jump matrices that are uniformly close to 
the identity matrix. The RH problem for R is also normalized at infinity. 

The proofs of Theorems 12.71 and 12.81 are in the final Section 9. 

3. The first transformation Y h^ X 

The starting point is the RH problem ()1.26p with functions (ll.27p . 
In this section we introduce the first transformation Y ^^ X which makes 
use of the special form of the weight functions Wj^n- 

3.1. The main idea. The main idea behind the first transformation Y 1-^ 
X can be found in an unpublished manuscript of Bertola, Harnad and Its 
|16j . The starting point is the observation that wo^n{x)e^^^' satisfies a 
scaled version of the so-called Pearcey differential equation 

(3.1) p"'iz)=zpiz). 

Special solutions to the equation (j3.ip are given by the Pearcey integrals 

(3.2) p,iz)= f e-^'/4+-ds, J = 0,...,5, 

Jr, 

where the contours Tj are 

Fq = (—00,00), Fi = (ioo,0] U [0,00), 

(3.3) F2 = (ioo,0]U [0,-00), F3 = (-ioo,0]U [0,-00), 
F4 = (— ioo, 0] U [0, 00), F5 = (— ioo, ioo). 
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Figure 3.1. Contours T,; in the definition of the Pearcey intergals. 



or homotopic deformations such as the ones shown Figure 13.11 Each Tj 
is equipped with an orientation as shown in Figure 13. 1[ Since the Pearcey 
equation (j3.ip is a linear third order equation, there must be a hnear relation 
between any four solutions. Indeed, from the integral representations (j3.2p 
one can find for example, 



(3.4) 



Pniz) =Pi{z) -pi{z). 



The functions Wj^n from (jl.27p are expressed in terms of the Pearcey 
integral po a-nd its derivatives as follows 



(3.5) 
(3.6) 
(3.7) 



W2,n{z) 



n-i/2e-"^(^)p'o(n3/Vz), 
n-3/4e-"^(^)p(;(n3/Vz). 



Now we can illustrate the main idea behind the transformation Y 
Consider the matrix-valued function X defined by 



X. 



(3.8) 
where 

(3.9) 



Xiz) = Yiz) 







DnP-\n'^/^Tz] 



Dr, 




n 



1/2 






p 



'Po Pj 
P'o P'j 

kpo p'; 



and j and k are such that po, pj, and pk are linearly independent solutions of 
the Pearcey equation, so that P is indeed invertible. To compute the jumps 
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for X we note that (I3.5I1-(I3.7I) can be written as 



'1\ (wo,n{zy 

(3.10) D-'P{n'/\z) I = e"^(^) u;i,„(z) 

,0/ \W2,n{z), 



SO that for x G M, by (fL26|l and pIHll . 
(3.11) 



x-{xy^x+{x) 



1 



/I 






(wo,n(a;) tt'i,n(a;) 'W2,n(a;)) D„P ^{n'^/^Tx] 
h 



g-ny{x) o\ 

1 



1 
1/ 



The jump matrix has only one non-trivial upper triangular 2x2 block. This 
2x2 block has a varying exponential in its 12 entry and so is precisely of 
the form that appears in the RH problem for orthogonal polynomials on the 
real line with respect to varying exponential weight [3011331 ITT] . So the jump 
matrix for X has a form that seems very promising for asymptotic analysis. 
The complication however, is that the asymptotic behavior of X becomes 
rather involved due to the nontrivial asymptotic behavior of the Pearcey 
integrals Pj{z) as z — > oo. We also have to deal with the Stokes phenomenon 
which implies that the asymptotic behavior may be different as z — > cxd in 
different sectors in the complex plane. For example, by classical steepest 
descent methods one can prove as in |20] that 



(3.12) po 

oo with Rez > 0, and 



(^) = y^z-i/3eK/^(l + 0(z-4/3)) 



as z 



(3.13) 



Po{z) 




i/SgK-^) 



4/3 



(l + 0( 



-4/3 



)) 



as 2 ^ OO with Rez < 0. Here and throughout the paper we use principal 
branches of the fractional powers, that is the branch cut is along the negative 
real axis. So po has a different asymptotic behavior as z — > oo in the right 
half-plane and in the left half-plane. Also the asymptotic formulas (|3.12p - 
()3.13p are only uniformly valid if we stay away from the imaginary axis. 

To deal with the different asymptotic behaviors we define P differently 
in each of the quadrants and define the transformation Y ^^ X accordingly. 
This procedure will introduce a new constant jump on the imaginary axis, 
while still simplifying the jump matrix on the real axis. Also the asymptotic 
condition in the RH problem for X takes a nice form. 
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3.2. A RH problem for Pearcey integrals. We denote the four quad- 
rants in the complex plane by /, II, III and IV , respectively. We define P 
in the different quadrants by 



(3.14) P 



(3.15) P 



Po 


-P2 


-P5 


P'o 


-P'2 
-P'i 


-P'i 


Vo 


-Pi 


-P5 


p'o 
/o 


-P4 
-P'i 


-p'i 



in /, 



P 



in ///, P 



Po 


-Pl 


-P5 


p'o 
p'o 


-p'i 
-p'i 


-P', 
-p'i 


'po 


-P3 


-P5 


p'o 
p'o 


-Ps 

-p'i 


-P5 
-Pi 



in//, 



mlV. 



Hence each column of P contains a particular Pearcey integral pj , see (13. 2p , 
and its derivatives. We use po in the first column and p^ in the last column. 
In the middle column we use different Pearcey integrals in the different 
quadrants. 

Then P is analytic in C \ (M U iM) with a jump on M U iM. By linear 
relations such as ()3.4[) one can easily obtain jump conditions for P. In the 
following lemma we show the RH problem that is satisfied by this function 
P. 

Lemma 3.1. We have that P satisfies the following RH problem: 



(3.16) 



P is analytic in 

P+{z)=P^{z) 

P+{z)=P.{z) 



\ 



MUil 
0^ 



1 

-1 
-1 

1 







1 

0^ 


1, 



z e 



zG ill 



Piz) 



_ /z-V3 

27r(/ + 0(z-2/3)) 1 

V 




BjQ.j{z), 



as z —> cxD in the jth quadrant, 



where the constant matrices Bj are given by (we use lo 



(3.17) Bi 



(3.18) ^3 




Bo 



Ba 




1, 



and the diagonal matrices Qj{z) for j 
(3.19) 

Gi(z) = I e^i(^) , Q2{z) 
e^2(^)/ 



, 4 are given by 



i{^) 




e^3(^) 











Qd2{z) 











e^iC-) 
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(3.20) 

/e^2(2) \ 

G3(z) = e^3W , Q^{z) -- 
V e^i(^)/ 

with functions 6j defined as 

(3.21) ej{z) = -uj^z'^/^, i = 1,2,3, -7r<argz<7r. 

The asymptotics for P are uniform as z ^ oo in any region provided we 
stay away from the axis, that is, it is uniform as z —> oo in the region 

(3.22) dist ('arg z, ^z) > e, 

for some e > 0. 

Proof. In [20] a detailed discussion of a RH problem for Pearcey integrals 
can be found. The lemma is proved by following these arguments. However, 
we do wish to make a few remarks on the asymptotics of P. From classical 
steepest descent arguments one can prove as in [20] that 

(3.23) p,{z) = Mz~'/hh'^\l + 0{z-^/')) 

V 'J 

(3.24) p2{z) = a;^/^z-i/3e!-^/='(i + 0{z-'/')) 

V 'J 

(3.25) p,{z) = u.^z-y\i-'^"\l + 0(z-^/3)) 

as z — > oo in the first quadrant. In [20| these formulas (with a different 
numbering of the functions pj) are stated with an error term ©(z^^'^). 
However, our case is a special case of the more general case considered in 
|20j . By analyzing the proof for our special case we obtain that the error term 
is indeed of order 0{z~^''^). From these asymptotics (and the corresponding 
ones for the derivatives of the Pearcey integrals) we obtain by (13.14p that 

pr- /z-i/3 \ 

(3.26) P{z) = J-\ 1 (l + 0(z-4/3))5iei(z), 

^ ^ V z^l^j 

as 2 ^ oo in the first quadrant. Since 

Vi/3 \ 

(3.27) I 1 {l + 0{z-^'^)) 

z'^/'y 

'z-^l^ \ 
(1 + 0(^-2/3^)1 1 , asz^oo, 

zi/3/ 
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we obtain the asymptotic behavior of P in the first quadrant. The other 
quadrants can be dealt with in a similar way. 

Finally, because of the Stokes phenomenon, the asymptotics for pQ are 
only uniform if we stay away from the imaginary axis. Similarly, the asymp- 
totics for p5 are uniform if we stay away from the real axis. Hence the 
asymptotic behavior is uniform if we stay away from both the real and 
imaginary axis and let z — > cxd so that it remains in the region given by 
^^?F2\) for some e > 0. D 

3.3. The inverse transpose of P. In the transformation y i-^ X we will 
use 



(3.28) 



Q = P 



-t 



Then Q satisfies a RH problem which is easily obtained from the RH problem 
satisfied by P. 

Lemma 3.2. Define Q = P^*. Then Q satisfies the following RH problem: 

Q is analytic in C \ (M U iM) , 

'I 
Q+iz) = Q.{z) I 1 1 

1 

1 
Q+iz) = Q^{z) I 1 1 

^0 1^ 

/zi/3 

^(i + o{z-y^)) [ 1 \A,e,{z)-\ 

V z~^/^^ 
as z ^ oo in the jth quadrant, 



(3.29) 



z G 



z G 



Qiz) 



where Aj = B- , and the Bj and Qj, j = 1,2,3,4, are given by ()3.17p - 
(13:20]) . 

The asymptotics for Q are uniform as z — > oo in any quadrant such that 
([3:22]) holds for some e > 0. 



Proof. This follows immediately from (j3.28p and Lemma |3. 11 

The constant matrices Aj are given explicitly by (with lo = e^'^''^). 



D 



(3.30) 



(3.31) 



A,=^ 



A. 




U 1 UJ^ 

Ao = ^{ 1 1 1 



V3 



' ijj^ 1 LO 
1 1 1 



Aa 



V3 



V3 



vW"' 1 LO 

-1 U? 
-1 1 



UJ 1 LO j ^^ \— 1 LO UJ 

The prefactor i/\/3 is so that all matrices Aj have determinant 1. 
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3.4. The transformation Y h^ X. For clarity we define the transforma- 
tion Y ^^ X m. two steps. First we define Y ^^ X and then X ^^ X. 
We define X by 



(3.32) 



X{z) 



Cnr^^'U) DnQ(rv>I^Tz\ 



where -D„ and Q are given in (j3.9p . (|3.28p . and the constant prefactor C„ is 
given by 

^r-i/3 
(3.33) C„ = -y/^ I 1 

ri/3 



We also define 
(3.34) 



X{z)=X{z) 



1 







ej(n3/Vz) 

for z in the jth quadrant, where Gj is given by (J3.19p - (l3.20p . 

The 4x4 matrix valued functions X and X are analytic in each of the 
four quadrants with jumps across the real and imaginary axes. 

Lemma 3.3. We have that X satisfies the following RH problem: 
(3.35) 

X is analytic in C \ (M U iM), 

fl e-"^(^) 0\ 
1 

11 

1/ 

1 0\ 
10 



X+(z)=X_(z) 



z G 



X.(z) 



X-{z) 



110 

\o ly 



zG ill 



X(z) = (/+0(z-2/3)) 



n+1 

z — 3- 















n 

Z 3 





\ 





1 

AjQ^^{n^l^Tz) 





n-l 
Z 3 / 

as z ^ oo in the jth quadrant 



The asymptotics for X are uniform, as z ^ co in any quadrant provided we 
stay away from the axis, that is, it holds uniformly as z —> cxd with (I3.22p 
for some e > 0. 



Proof. Let z G M. Then we obtain from (fOHl) and (f3^32]) that 



(3.36) XZ 



l( \Y ( \-f^ {W0,ni^) Wi^n{z) W2,n{z)) DnQ+ 

[z)Ji+(z) - y^^ Q_{n^/^Tz)-^Q+{ny^Tz) 



{n^/^Tz) 
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Since Q = -P~*, we can use (|3.1ip to find 

{wQ,n{z) Wi,n{z) W2,n{z)) D r,Q +{1^^ ' ^ T z) = (e"" ^(^) O) . 
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Using also the jump condition on the real axis in the RH problem ()1.26p 
for Q, we obtain the jump matrix for X on the real line as given in the RH 
problem (j3.35p . 

The proof of the jump on the imaginary axis is similar (even simpler). 

Finally we check the asymptotic formula for X. From the asymptotic 
behavior of Q given in (j3.29p and the definitions (|3.9p and ()3.33p . of D„ and 
Cn, we obtain as z ^ oo in the jth quadrant, 



(3.37) 

(z^l'^ 

DnQirv'/^Tz) = C-\I + 0{z-^/^))\ 1 \ AjQ-^rv'/^Tz). 

\ z-^/'\ 



From the asymptotic condition in the RH problem (jl.26p for Y and the 
definition ()3.32p . it follows that 



(3.38) X{z)={\ ^\{l + 0{z-^)) 



(z^ 

z-"/3 



Vo 



\ 



^-n/3 

z-"/V 



1 







DnQin^l^TZ 



Combining this with (j3.37p gives the asymptotic condition for X in (j3.35p . 

D 



The transformation (j3.34p from X to X has the effect of simplifying the 
asymptotic condition in the RH problem, but the jumps on the real and 
imaginary axis are more complicated. 

Lemma 3.4. We have that X satisfies the following RH problem: 



X is analytic in C \ (M U iM), 
X+{z) = X_{z)Jxiz), zGMUi: 

/z" 



(3.39) < 



X(z) = (/ + 0(z-2/3)) 



n + l 

z — 3- 






Tl 

Z 3 



\ 




z~3 

\0 z'^J 

as z ^ oo in the jth quadrant 



1 

A, 
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where the jump matrices Jx are given by 
(3.40) 





Jx{z) 
for 2 € IE 
(3.41) 





Vo 

-, and 



Jx{z) 



1 



/I 





Vo 













\ 



4/3 



^:fn3^ir4/3|^|4/3 





„3V3i^4/3|^|4/3 


0\ 




1 


g±n3firV3|^|4/3 ^ 





ij 



for z e iR± . 

The asymptotics for X are uniform as z ^^ oo in any quadrant provided 
we stay away from the axis, that is, it holds uniformly as z ^ oo with (j3.22p 
for some e > 0. 

Proof. The asymptotic condition is clear from (j3.34p and the asymptotic 
condition in the RH problem (j3.35p for X. 

Let z G IR+. Then it follows from (j3.34p . the jump condition in the RH 
problem (j3.35p . and the expressions for 0i and 64 from (J3.19p - (j3.20p that 



(3.42) 

XZ\z)X+iz) 



1 

e^\n^/*z) 



(l e-'"^(^) 0\ 

1 

11 

1/ 



1 



\0 







^ex{n^/^Tz)-e2{n^'^Tz) 





1 

Qi{n^l^z] 



\ 



1 

^ei(n='/'4r2)-e2{n3/4r2) 



Inserting Oj{z) 



|c^j'z4/3 



we arrive at the jump matrix (j3.40p for z £ 



The proof of (j3.40p for z £ R_ is similar, taking into account that the 
functions 6j have their branch cuts along the negative real axis. 

The proof of the jump matrix (j3.4ip for z G iM is similar. D 

The fact that the asymptotics for X is not uniform up to the axis is a 
complication that will be automatically resolved during the steepest descent 
analysis. 

The RH problem for X reveals to some extent how the equilibrium prob- 
lem comes into play. We see that the jump matrices in the RH problem 
locally have the structure of jump matrices in a 2 x 2 matrix valued RH 
problem. The measure //i will be used to handle the upper left 2x2 block 
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in the jump matrix (I3.40p on M and the measm'e /U3 will handle the lower 
right block. The measure ^2 is supported on the imaginary axis and will 
handle the middle 2x2 block in the jump matrix (I3.41[) on the imaginary 
axis. 

Thus /xi will act on the first and second rows and columns, ^2 and the 
second and third, and /i3 and the third and fourth. Then we also see that 
/ii and ^2 interact in some way, since they both act on the second row and 
column. Similarly, 1x2 and /X3 will interact. On the other hand, the measures 
//I and /i3 do not interact with each other, since they act on different rows 
and columns. This is reflected in the equilibrium problem where we have 
that the mutual logarithmic energies I(^i, /X2) and I(/i2, /U3) are present, but 
not the mutual logarithmic energy of ^1 and /i3. 

The potential V{z) in ()3.40p comes with an extra term — |r^'^|z|'*'^. To- 
gether they appear in the equilibrium problem as the external field on /ii. 
The lower 2x2 block in (l3^ and the middle 2x2 block in (f3^ have 
rapidly oscillating entries on the diagonal. This is reminiscent of what hap- 
pens in the steepest descent analysis of the RH problem for orthogonal 
polynomials after the application of the gf- functions |33j . 

We will explain at a later stage why the restriction 122 ^ cr plays a role. 
Here we can only say that it is partly explained by the triangular structure 
of the jump matrices (13.4ip and ()3.40p . Indeed the non-trivial 2x2 block 
in (j3.4ip is lower triangular, while the 2x2 blocks in (j3.40p are upper 
triangular. 



4. Analysis of the equilibrium problem 

To prepare for the next transformation of the RH problem we need the 
equilibrium problem for the energy functional Ey- In this section we prove 
Theorem 12.31 which gives existence and properties of the equilibrium mea- 
sures Hi, 1^2 and /i3. We also prove Propositions 12.41 and 12.61 A basic tool 
is the balayage (sweeping out) of a measure. 

4.1. Balayage. The balayage of a finite measure v onto a closed set K with 
positive capacity is a positive measure v on K with \\v\\ = \\V\\ and with the 
property that there exists a constant C such that 

(4.1) V'ix) = U^{x) + C, q.e. x e K, 

where q.e. means quasi-everywhere, that is, with the exception of a set of 
zero capacity. We will be interested in cases where K is either the real or 
imaginary axis, or K consists of two unbounded intervals along the imagi- 
nary axis 

(4.2) Kc = (-ioo,-ic] U [ic,ioo), c> 0. 
In these cases we have C = for the constant in (14.111. 
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We use the notation u for the balayage oi u ii K is understood. To 
emphasize K we write 

(4.3) d = Bal{u,K). 

Each balayage measure can be written as the integral over balayage mea- 
sures of Dirac-delta measures as follows, 

(4.4) Bal(i/,i^) = I B&\{5^,K) dv{z). 

Hence it is sufficient (for our purposes) to calculate Bal((5z,i^) for several z 
and K. 

Example 4.1. If i^ = M is the real line, and y ^ 0, then the balayage 
measure of 5iy onto M has the Cauchy density 

dBalW,.R)(.) ^ l\y\_ ^, ^ ^^ 

dx vr j;^ + y^ 

For i^ = iM and x 7^ we have similarly 
1^ ^, dBay jR) ^ !_ 1^^ ^ ^ ^^ 

|dz| vr |z|^ + x^ 

These results are well-known. They are also the limiting case c -^ 0+ of 
3]) in the next example. 



Example 4.2. Let Kc = (— ioo,— ic] U [ic, ioo) with c > 0. If y G {—c,c) 
then the balayage of 5iy onto iCc has the density 



2 



(^ ^, dBal(i. «,) ^ 1 yg^ ^ , ^. 



|dz| 7r|z-iy|/| 



zp — c^ ' 



If X € M then the balayage of 6x onto -fCc has the density 

dBal(^.,i^,) _ 1 |z|^/?+^ ^^^ 

\dz\ ^(|z|2 + x2)v/|z|2-c2' 

In all cases we have that the constant C in (j4.ip vanishes. 

These expressions can be proved by simple contour integration. We omit 
the proofs. It will be important for us that the densities (j4.7p and (j4.8p are 
decreasing as \z\ increases, with a decay rate 0(|2;|~^) as \z\ -^ 00. 

In the following we will see how the notion of balayage is related to the 
equilibrium problem. 

4.2. Equilibrium problem for 1^3. If we fix two of the measures Uj, j = 
1, 2, 3 we can consider the equilibrium problem with respect to the remaining 
measure only. If 1^1 and 1^2 are given, then the equilibrium problem for 1/3 is 
to minimize 



(4.9) I{u) - I 



W^du 
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with respect to all measures z^ on M with total mass 1/3. This problem 
only depends on 1/2- Assuming that 2^2 has finite logarithmic energy and 
total mass 2/3 (as in conditions (a) and (c) of the equilibrium problem), the 
minimizer for ()4.9p is given by z^s = 2 Bal(t'2,]K) since it satisfies W^^ = 2^^^^ 



on M, which is the variational condition for (j4.9p . and it has the correct total 
mass 1/3. So we have 

Lemma 4.3. Suppose vi and V2 are fixed so that conditions (a), (b), and (c) 
of the equilibrium problem are satisfied. Then the measure v^ that minimizes 
EY{vi,V2,u-i) subject to the condition (d) is given by 

(4.10) z.3 = ^Bal(i/2,M). 

Thus by ()4.5p and ()4.4p . we see that 1/3 has a density given by 

——{x) = - -^ : 

ax TT J x"^ + \z\ 



"(^) = - / ^2 I U|2 d'^2(^)- 



4.3. Equilibrium problem for i^i. For given 1/2 and z/3 having finite log- 
arithmic energy, the equilibrium problem for z/i is to minimize 



U^'^ix) dz^(x) 



(4.11) I{u) + / (v{x) - ^r4/3|x|4/3 _ U"' (x)] 

among all probability measures v on M. This is a usual equilibrium problem 
on M with external field 

y(a;)-^r^/3|a,|4/3_^^2(^). 
4 
The term — C/^^ jj^ i^j^g external field attracts the measure 1^ towards 0. 

Lemma 4.4. Suppose 1^2 o-nd 1^3 are fixed so that conditions (a), (c), (d) 
and (e) of the equilibrium problem are satisfied. Then the measure vi that 
minimizes £'y(i/i, 1/2, t's) subject to the condition (b) is the minimizer for 

Km . 

It has the following properties: 

(a) The support of ui is contained in the convex hull of the support of the 
equilibrium measure on M. in the external field V{x) — jT^'^\x\^'^ . 

(b) If 1^2 is so that there exists a constant c > such that S{a — V2) C K^ 
where Kc is given by (14. 2[) . then 

4 ' ' 
is real analytic on M. 

Proof, (a) We use the fact that the support of the equilibrium measure 
with external field Q on R is the closure of those points where weighted 
polynomials e~"''^P^, degP„ < n, n G N, take their maximum modulus [7T| 
Chapter III.2]. 

Assume the equilibrium measure in the external field Q{x) = V{x) — 
|r'^'3|j;|4/3 ^g supported in [— X, X] for some X > 0. Hence the maximum 
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modulus of e""'^^^, degP„ < n, is attained only in [— X, X]. Since U^'^{x) = 
J log ,_, (li'2{z) is even on M and strictly decreasing as \x\ increases, it 

then follows that the maximum modulus of e^^"^ e~^^ P^ = e~"'^'^~'^'"^^Pn, 
deg Pn ^ n is attained on [—X, X] only, which proves part (a) of the lemma, 
(b) Suppose 1^2 and c > are such that S{a — V2) d Kc- Then 

(4.12) U'"'{x) = - f log\x-z\da{z)- [ \og\x - z\di^2iz), 

and the second term is real analytic for x G M. For the first term in (j4.12|) 
we have by the definition (j2.7p of a, 

(4.13) - I" log \x - z\da{z) = -^r^/' y"'log(x2 + y^)y'/^dy 
whose derivative with respect to x is 

(4.U) -f.-f;|^*. 

71" Jo x^ + y'^ vr Jc x^ + V 

The second term on the right-hand side of (j4.14p is real analytic on M. In 
the first term we introduce the change of variables y = \x\s to obtain 

Because of the standard integral 

(4-15) / T— ^ds = - — , -l<p<l, 

Jo 1 + s^ 2cos(p7r/2) 

which for p = 1/3 gives the value 3 yovr, it then follows that 

TT Jo X'^ + y^ 

Combining this with (|4.12p . (|4.13p . and (|4.14p . we obtain 

-—U'''^{x) = -r^/^x|x|~^/^ + "real analytic function", 
dx 

which after integration completes the proof of part (b) of the lemma. D 

Lemmas 14.31 and 14.41 show how the solutions of the equilibrium problems 
for ui and u^ are determined by the measure f2- Next we study the converse: 
given ui and ^3, what are the properties of the solution of the equilibrium 
problem for 1^2? 
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4.4. Equilibrium problem for U2. Now suppose that we are given vi and 
1/3 satisfying the conditions (a), (b), (d) of the equihbrium problem and we 
wish to minimize Ev{vi, 1^2, '^3) with respect to z^2- The equilibrium problem 
for 1^2 is to ininimize 

I{u)- f {U"' + U"^') du 

with respect to all measures v on iM with total mass 2/3 and satisfying the 
constraint v < a. 

If the constraint v < a were not present, then simple balayage arguments 
would give that the minimizer is given by 

-Bal(z.i + z.3,»). 

Note that this measure has indeed the correct total mass 2/3. However, this 
measure will violate the constraint. Indeed, the density of the balayage of 
any measure on M onto iM has its maximum at 0, and is strictly decreasing 
for 2; G iM if |z| increases, see also ()4.6p . Since the density of a (|2.7p vanishes 
at 2; = and is strictly increasing if \z\ increases, it is clear that the density 
of the balayage measure lies strictly above the density of a on a non-empty 
symmetric interval around on iM. 

Lemma 4.5. Suppose that vi and 1/3 are fixed so that conditions (a), (b), 

and (d) of the equilibrium problem are satisfied. Then the measure 1^2 that 

minimizes Ev{i'i,i'2-,i^3) subject to the conditions (c) and (e) exists, and 
there is a constant c > such that 

supp((T — z^2) = Kc := (— ioo, — ic] U [ic, ioo). 

Moreover, V2 satisfies 

iu'^^{z) = U-^{z) + U''Hz), zGiM\(-ic,ic) 
^' ^ [W^iz) <U''^{z) + U''^i{z), ze{-ic,ic). 

and a — U2 vanishes as a square root at ibic. 

Proof. We use the iterated balayage algorithm of [53j . Write 

U2fi ■■= -Bal(z/i + z/3,iM) . 
As already noted above, there exists cq such that 

'S'((z^2,o -cr)^) = [-ico,ico] 

where {1^2,0 — (^)~^ denotes the positive part of the signed measure 1/2,0 ~ cr- 
It then follows from the saturation principle of [39] that the minimizer V2 
(if it exists) is equal to the constraint a on [— ico,ico]. 
Then we define 

(4.17) 1/2,1 := min((T, z/2,0) + Bal ((1^2,0 - o-)+, Kco) • 

In other words, we balayage the part of f2,o that lies above a onto the part 
of the imaginary axis Kcq where the constraint is not (yet) active. Then 1^2,1 
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has a density for z G K^^ that is strictly decreasing as |^| increases, see also 
(|4.7p . Then there exists c\ > cq so that 



5((i^2,i - 0")"^) = [-ici, -ico] U [ico,ici] 
and we define 

(4.18) 1^2,2 := min((T, z^2,i) + Bal ((1^2,1 - o-)+, Kci) • 

Continuing this way we find an increasing sequence (c^) and a sequence 
{t^2,k) of measures of total mass 2/3. Clearly, the Ck must be such that 

(4.19) (T{{-ick,ick]) < -. 

Hence the sequence converges and we denote the limit with c. It follows as 
in [53] that (2^2, fc) converges weakly to a measure f2- In fact, after each step 
in the iterated balayage process we have that 

(u-^-'^+^{x) = U-^-'^{z), zGiR\(-icfc,icfc), 
^ ■ ' [W'^^^+^ix) < U'^^^'^iz), z G (-icfc,icfc). 

From these equations and the fact that 2C/''2,o (2) = [/^i (^z) + U"^ (z) for 
all z £ iM, we can easily deduce (j4.16p . Hence 1^2 satisfies the variational 
conditions (with strict inequality) that uniquely characterize the minimizer. 
From the above arguments it also follows that the density of 1^2 for z G Kc 
is strictly decreasing as |2;| increases. Therefore the density of z^2 hes strictly 
below the density of a on (— ioo, — ic) U (ic, ioo) and there is a square root 
vanishing of the density of a — 1^2 at itic. D 

4.5. Proof of uniqueness of the minimizer. We start by rewriting the 
energy functional as 

(4.21) Ev{i^i,U2, 1^3) = ^I{i^i) + ^/(2z^i - 31/2) + jl{i^2 - 21^3) 



+ f fv{x)-^T''/^\x\'AdMx). 



If two measures ly and fi have finite logarithmic energy and j di^ = J dfi, 
then 

(4.22) I{u - /u) > 0, 

with equality if and only if z/ = fi. This is a well-known result if vi and 1^2 
have compact supports [7l]. For measures with unbounded support, this is 
a more recent result of Simeonov [72], who obtained this from a very elegant 
integral representation for /(z^ — fi). 

It follows from (I4.22P that the energy functional ()4.2ip is strictly convex. 
The minimizer is therefore uniquely characterized by the Euler-Lagrange 
variational conditions listed in Proposition 12.41 
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4.6. Proof of the existence of the minimizer. We rewrite Ey again as 
in ([421]) . From (14221) it follows that 

(4.23) I{2ui - 3z^2) > and I{u2 - 2^3) > 0. 

Combining this with ()4.2ip leads to 



dz^i(x) 



(4.24) Ev{ui,iy2,i^3) > \l{yi) + j {v{x) - ^r^/^lxl^/^) 

and this implies that Ey is bounded from below. 

A natural approach to prove the existence of the minimizer is now the 
following. Since Ey is bounded from below we can approach the infimum by 
a sequence of measures satisfying the conditions of the equilibrium problem. 
If we can prove that this sequence is tight, then we have that the sequence 
converges weakly to the minimizer of the Ey. However, it seems complicated 
to follow this procedure. The main difficulty is the fact that the measures 
f2 and f3 will have unbounded supports with rather fat tails. It turns out 
that their densities behave like 0(|z|~^'^) as \z\ -^ 00. 

Our method of proof consists of restricting the measure f 2 by an additional 
constraint 

(4.25) V2 < K(Tp 
where (jp is the measure on iM with density 

(4.26) dap{z) = -—\dz\, z G iM, 

with p G (1, 5/3) and K > Q is a. suitable constant that will be determined 
later on. So we consider the equilibrium problem with extra constraint (j4.25p 
and we prove that for the minimizer the extra constraint is not active. 
The special role of the value p = 5/3 is made clear in the following lemma. 

Lemma 4.6. Let 1 < p <2. If 1^2 < KcTp and u^ = 2 Bal(z^2)M) then 

with constant 

(4.28) C„ = . --. 

^ ^ ^ 2sin(p7r/2) 

We have Cp < 1 if and only if p < 5/3. 
Proof. By direct calculation 

dus If \z\ . I .^ ^ r \y\ ^ A 



dx 2Tr J x^ + |zp 2tt J_^ x"^ + y"^ \y\P 

= — / -o- — 2^y, X G 
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Make the change of variables y = \x\s and use ()4.15p to obtain 

di/3 K f^ s^-P , K vr 

< , , / Kds 



dx Tr\x\P Jq 1 + s^ 7r|a;|P 2cos((1 — p)7r/2) 

K 1 

X G M. 



|x|P2sin(p7r/2)' 
This proves the lemma. Note that indeed Cp < 1 if and only if p < 5/3. D 

We may iterate the above argument and we obtain 

Corollary 4.7. If 1^2 < Kap and v^ = ^Bal(i/2,M), then 

(4.29) - Bal(zv3, iM) < ClKap. 

We fix p € (1, 5/3) and we show first what constant K > Q to take. Let 
X > be such that the support of the equilibrium measure in the external 
field V{x) - |r^/3|^|4/3 jg contained in [-X,X]. 

Prom the explicit formulas (j4.7p and (j4.8p for the balayage of point masses 
in [— ic, ic] U M onto K^, we can easily find the bounds 

dB&\{5,y,K,) ^ 1 c 

m ~Ti2 — 2' ^ ^ -'^c, 

\az\ TT \z\'^ — c^ 

and 

dBal(J^,i^e) ^ 1 ^/?T? 

U I — I 12 2 ' ^ ^ f^" 

|dz| TT |zp — C^ 

From this it follows that 

dB a\{p,Kc) ^ 1 Vc2 + X 2 

|dz| ~ vr I^P — c^ 

for every probability measure p on [—X, X] U [— ic, ic]. 

Let Co be such that cr([— ico,ico]) = 2/3. It follows from (j2.7p and (j4.30p 
that there exists a constant L = Lx so that 

/.O.N . / d^ dBal(/9,/s:cn)\ L 

(4.31) mm -— , rrv^ ) ^ V^^ ^ ^ if , 



(4.30) ""T; "^ <Z ^r 2 - -e^- 



dzl ' Idzl / \z\P 



'CO 



for every probability measure p on [— X, X] U [— ico,ico]. Since Cp < 1 we 
can then find a constant K large enough so that 

(4.32) C^K + \l<K. 

This is how we choose K and this determines the extra constraint Ka^. 

We now claim that there exists a vector of measures (i/J', i^^j '^s) that 
minimizes Ev{i^i, ^2) i^s) subject to the conditions (a) — (e) of the equilibrium 
problem of Definition 12.21 and subject to the additional constraint 

(4.33) 172 < Kap 
with K satisfying (143211 . 
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We already found that Ey is bounded from below. Therefore there exists 
a sequence {1^1,71,1^2,^1^3,71) such that 

(4.34) Evil^l,n, V2,n, l^3,n) < - + ml Ey 

where both the sequence and the infininiuni are taken with respect to condi- 
tions (a)-(e) of the equilibrium problem and (I4.33p . By Lemma 14.31 we may 
as well assume that 

(4.35) I/3,n = ^Bal(l/2,n,M). 

Moreover, we take i^i^„ as the minimizer of E{i'i, 1^2,71,1^3,71) subject to con- 
dition (b). By (j4.25p and Lemma 14.61 the sequences 1^2,™ and 1^3,71 are 
tight. Moreover, by Lemma 14.41 and the choice of X we have that the mea- 
sures 1^1,71 are supported in [— X, X]. Therefore there exists a subsequence 
{'^i,kn, ^2,k„, ^3,fc„) that converges weakly to a vector of measures {ul, i^l) ^3)- 
This is the minimizer that we seek for and hence we proved the claim. 

Our final goal is to show that the extra upper constraint Kap is not active. 

Given the measures v^ and 1^3 from the minimizer of the equilibrium 
problem with extra constraint, we denote by 1^2* the measure that minimizes 
the energy functional -By (i^*, 1^2, 1^3) subject to the constraint 1^2 < cr only. 
By Lemma [4. 5 1 we then have that 1^2* is equal to a on some interval [— ic, ic]. 
We clearly have c < cq. It also follows from properties of balayage that 

4*\K^^<^Bal{,.l + i^*3,K,,) 

1 
2 
where 



(4.36) < - Bal(i^3* , iM) + Bal(y9, K,, ) 



p = iz.i* + iBal(z.3*,i]R)|j_,,,,.,,,] 

is a measure on [— X, X] U [— ico,ico] with total mass < 2/3. Recalling also 
that 

1^2 <o- 
we find by Coroharyllll K3T\\ . and (f06]) that 

diyi* ^oK 2 L 

TTT<Cl-— + --— zGiM. 

|dz| - P\z\P 2,\z\P 

Then it follows from ()i:32|) that 

(4.37) ?^ < r^, ^ e iM. 

Thus V2* would also be the minimizer if we impose the extra constraint 
V2 < KcTp. It follows that 

V2 = V2 

and since the inequality in (j4.37p is strict, it follows that z^| does not feel 
the extra constraint. 
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Thus (z^j*, 1^2 ^J^s) is ^Iso t^^ niininiizer of the equihbrium problem of Defi- 
nition [22] (without extra constraint), and existence of the minimizer is now 
proved. 

4.7. Proof of Theorem 12.31 and proof of Proposition 12.41 

Proof. We have shown that a unique minimizer (/xi,;U2,/i3) for the equihb- 
rium problem exists. 

Proposition 12.41 and parts (a)-(c) of Theorem 12.31 follow by taking the 
minimizer (/xi,/i2, Ais) for the vector {I'l, 1^2,1^3) in Lemmas 14.51 14.3] and 14.41 

Part (d) of Theorem 12.31 follows from the symmetry in the equilibrium 
problem (note that V is even) and the uniqueness of the minimizer. D 

4.8. Proof of Proposition 12.61 

Proof. The measure //i is the unique minimizer of the energy functional given 
in ()4.1ip . Since the external field is symmetric, it follows from [71 1 Theorem 
IV. 1.10 (f)] that d^i(t) = dJl{t'^)/2 where Jl is the unique minimizer of the 
energy functional 

(4.38) 

// log ]_ du{x)du{y) + 2 / (v{^/^) - t^I^x^I^ - V^^{^)\ di/(x), 

among all probability measures v with support on [0, 00). From 

(4.39) -U^^ {^) = A f log(x + \z\^)dfi2{z) 

it follows that x >-^ —U''''^{^/x) is a convex function on [0,oo). Also x 1— > 
—t'^'^x'^'^ is convex on [0, 00). Hence, since we assumed that x 1-^ V{^/x) is 
convex as well, the external field in the equilibrium problem for fl is convex. 
Thus the support of fi consists of one interval [a, b] with < a < 6, see 
e.g. |71t Theorem IV. 1.10]. Hence the support of fii is either one or two 
intervals, depending on whether a = or a > 0. D 

5. A RiEMANN SURFACE 

5.1. A four-sheeted Riemann surface. We use the solution of the equi- 
librium problem to construct a Riemann surface as follows. Start with the 
four sheets TZj defined as 

, . 7^l = c\5(^l), 7^2 = c\(5(^l)u5(a-/U2)), 

^^■') 7^3 = c\(5((J-^2)u5(M3)), 7^4 = c\5(/i3). 

The four sheets are connected as follows: T^i is connected to Tl2 via S{fii), 
7^2 is connected to 7^3 via S{a — ^2) and TZ^ is connected to TZ^ via S{^3), 
every connection is in the usual crosswise manner. See also Figure ISTTJ for a 
picture of the Riemann surface in the case where the support of fii consists 
of one interval. 

The Riemann surface is compactified by adding two points at infinity: one 
is on the first sheet, and the other point at infinity is common to the other 
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^3 



7^4 



Figure 5.1. The Riemann surface 7^ for the one-cut case 



three sheets. The Riemann surface has genus A^ — 1 if A^ is the number of 
intervals in the support of ni. 



5.2. The Cauchy transforms. The Cauchy transform F of a measure fj, 
is defined as 



(5.2) 



F{z) 



1 



z — X 



d/i(x), z £ C\ S{fi). 



Note that we slightly abuse the notion of Cauchy transform since 

ZVri J X — Z ZTTl 

is also called the Cauchy transform of fi. 

We are interested in particular in the Cauchy transforms of the measures 
/ii , //2 and /i3 that are the solution of the equilibrium problem of Definition 
12. 2[ We denote their Cauchy transforms by Fi, F2, F3, respectively. We 
use the Cauchy transforms Fj to construct a meromorphic function on the 
Riemann surface TZ. We recall the Sokhotski-Plemelj formulas according to 
which 

Fj,+ {z) + F,, _(z) = 2PV I -^d^i,{x), 

J Z X 

F,,+ (z)-F,,_(z) = -2^i^, 
where PV denotes the Cauchy principal value. 
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Lemma 5.1. The function ^ : \J-^TZj 

V'{z)-F,iz), 



C defined by 



(5.3) az) = I 



Fi(z)-F2(z)+r4/3zi/3^ 



F3{z) 



r4/3/ 



-z) 



1/3 



zGTZi, 

z £ IZ2, Rez > 0, 

z G 7^2, Rez < 0, 

z e 7^3, Rez > 0, 

z £ 7^3, Rez < 0, 

z £ 7^4, Imz > 0, 

z £ 7^4, Imz < 0, 



F2{z) 
F2(z)-F3(z)+r4/3zV3^ 

/las an extension to a nieromorphic function (also denoted by (,) onTZ. The 
nieromorphic function has a pole of order deg V — 1 at infinity on the first 
sheet, and a simple pole at the other point at infinity. 

Proof. The analyticity of £, on TZi and 7^4 is immediate. To obtain the 
analyticity of ^ on 7^2 and 7^3 we note that from the Sokhotski-Plemelj 
formula and from the fact that /X2 = o" on [— ic, ic] it follows that 



F2,+ {z) - F2,-{z) = -2m 



.d;U2(2) 



dz 



-27ri 



da(z) 
dz 



(5.4) 



^4/3^31^1^/3 ^-T^/3(zV3 + (_^)l/3) 



for z £ (— ic, ic). By the definition of ^ this implies that ^ has an analytic 
continuation to (— ic, ic) on both 7^2 and 7^3. Thus ^ is analytic on each 
individual sheet TZj . 

It remains to check the analyticity if we cross a cut and move from one 
sheet to another. That this holds, follows from the variational conditions 
for the equilibrium problem and the symmetry in the problem. We will only 
show that ^ is analytic when one crosses S(ni) since the other cases follow 
from similar arguments. 

Differentiating U^^ gives 

(5.5) 24-U'''ix) = -2PV f ^—dfiiis), x£S{^ll), 

ax J X — s 

which by the Sokhotski-Plemelj formulas can be written as 

d 



(5.6) 



-2— C/^i(x) 
ax 



Fi,+ (x)+Fi,_(2;), x£S{fii 



On the other hand, we can differentiate the right-hand side of the variational 
condition ()2.13p to obtain 

(5.7) 

d _ d 



2—U''Ux) = —U^'Hx) 
ax ax 



V' [x) + t\x\ / sgnx 



d_ 
dx 

'dx2 



/ log \x — s|d/i2(s) — V {x) + r|xp 3 sgn 
log(x2 - s2)d/i2(s) - V'{x) + t\x\^'^ 



sgnx 
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;-d/i2(s) — y'(x) + r|x| ' sgnx 



1 



-d^2(s) 



1 



X — S ' ' ' 'Z J X + s 
for X G S{fii). By symmetry of fi2 we obtain 



dfi2{s) — V' {x) + t\x\ ' sgnx 



- / d/i2(s) + - / — — d^2(s) 

2 J X — s 1 J X + s 



1 



-d)U2(s) 



Fofx) 



(5.8) 

x — s 

for X E M. Combining dSSD, (IS2D and dSSD leads to 

(5.9) Fi,+ (a;)+Fi„(x) = F2(x) + ^'(x) - t|x|^/^ sgnx, 

for X £ S{fii). This is exactly what is needed to prove that ^ is analytic 
when one crosses 5'(/Ui) and moves from TZi to 7^2- D 

We use .^j to denote the restriction of ^ to TZj . 

As a corollary to the lemma we find the following result about the as- 
ymptotic behavior of the Cauchy transforms. 

Corollary 5.2. There exists a constant a such that (with co = e 



2m 



") 



(5.10) 
(5.11) 

(5.12) 



Fi{z) = - + 0{z 



-2\ 



Foiz) 



Fsiz) 



l--^ + 0{z- 



2^ 

32 



5/' 



+ 0{z- 



Rez > 0, 

Rez <0, Imz > 0, 

Rez <0, Imz < 0, 

Imz > 0, 
Imz < 0, 



uniformly for z ^ oo. 

Proof. The statement about Fi is obvious, since /ii has a compact support. 
For z large enough define H as 



=(^)=6(^') 



vr 



vr 



- < arg z < 

D 



and then continued analytically into other sectors. Explicit calculations 
show that 



(5.13) 



■iz) 



Uz' 



esiz'), 



Uz'), 



— 7r/6 < argz < 7r/6, 
57r/6 < argz < vr, 
— vr < argz < — 57r/6, 
7r/6 < argz < 7r/3, 
— 7r/3 < argz < —n/Q, 
2vr/3 < argz < 57r/6, 
— 57r/6 < argz < — 27r/3, 
7r/3 < argz < 27r/3, 
— 27r/3 < argz < —tt/3. 
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Then H is analytic for |2:| > i? for R large enough with a Laurent expansion 

oo 

(5.14) E{z) = T^''^z + ^CjZ~K 

j=0 

Since ^j{—z) = —^j{z) we also have S(— z) = —E{z), and hence the even 
powers vanish in this expansion. 

Next we claim that we can reconstruct the ^j, j = 2,3,4 out of H again 
in the following way 

'h(zV3), Rez>0, 

, Rez < 0, Imz > 0, 

E{z'^/^Lo'^), Rez <0, Imz < 0, 

'h(zV3), Rez<0, 
H(z^/^a;), Rez > 0, Imz > 0, 
E{z'^/^uj'^), Rez > 0, Imz < 0, 

5(^1/3^2), lmz>0, 
E{z^^^uj), lmz<0. 

We leave the verification to the reader. Since we have that by (j5.3p 
(5.18) F2{z) = -Uz) - Uz) + V'iz) + r^/s^Vs 

for Rez > 0, we obtain 

1-C3 



(5.15) 

(5.16) 
(5.17) 



(5.19) 



Fo z) 



Cl 



+ 



C5 



+ 0(z''] 



for z — > 00 and Re z > 0. Since F2{z) = f j^d/i2(a;) and since ^2 has total 
mass 2/3, we see that ci = and C3 = 1/3. This proves the asymptotic 
behavior of F2 with a = C5 for Rez > 0. The asymptotics for the other 
sectors and other Cauchy transforms follow from similar arguments. D 

The solution to the equilibrium problem can be described in terms of the 
functions S,j as described in the following proposition. 

Proposition 5.3. With ^j = ^\fi-, we have that 
(5.20) dfii{z) = ^ (6,+ (z) - CiA^)) dz, 



27ri 
1 



Z £ 



(5.21) dfi2{z) = — (6,+ (^) - 6,- (2)) dz + da{z) 



zeil 



z G 



(5.22) d/i3(^) = ^ (6,+(^) - e3,-(^)) dz - l!^|z|V3dz, 

Proof. This follows by Lemma 15.11 and the Sokhotski-Plemelj formula 



(5.23) 

for z G S{fj,j). 



1 

'2^ 



(F,-+(z)-F,,_(z)) 



d_ 

di 



/^j( 



D 
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6. The second transformation X ^^U 

From now on, we assume that ^ui is one-cut regular. That is, ^i is sup- 
ported on one interval, and there are no singular points in the sense of 
Definition 12.51 We write 

with a > 0. 

6.1. Definition of the transformation. In the second transformation in 
the steepest descent analysis, we use the (^-functions associated with the 
measures /Xj that minimize the energy functional Ey- For j = 1 and j = 3, 
we define Qj as 



(6.1) gj{z) = j log(z - s) dnj{s), j = 1,3, 

with the principal branch of the logarithm 

log{z — s) = log Iz — s| + iarg(2; — s), arg(z — s) G (— v", n). 

We also define (72 by a similar integral, but since fi2 is supported on the 
imaginary axis, we make a different choice for the branch of the logarithm. 
We define (72 by 

(6.2) g2{z)= flog{z-s)dfi2{s). 

where we take the logarithm such that 

(6.3) log{z — s) = log\z — s\ + iai[g{z — s), aig{z — s) £ {—tt/2,3tt/2), 

for z G C \ iM and s G iM. Thus 52 is defined and analytic in C \ iM. 
The transformation X h^ U is now defined by 

(6.4) U{z) = L-^X{z)G{zTL'^, 
where G is given by 

(6.5) G(z) = diag(e-3i(^) e9i(z)-g2(z) f,g2(z)-g-i{z) ^93(2)) 

and L is given by 

(6.6) L = diag (e"^ 111), 

with a. the variational constant from (j2.13p . 

Since X and the Qj are analytic in C \ (M U iM) , it is clear that U is also 
analytic in C \ (M U iM). The jump matrices Ju in the RH problem for U 
are obtained from the jump matrices Jx by 

(6.7) Ju = L~"GI" JxG+i.", 

and we have to calculate them on the various parts of the real and imaginary 
axis. We will see that the jump matrices take a simple form that is suitable 
for further analysis. 
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Also the behavior at infinity simphfies under the transformation. How- 
ever, in contrast to other works (see e.g. [33j), the RH problem is not nor- 
malized at infinity after the use of the g-functions. So U{z) does not tend 
to the identity matrix as z — > oo, but instead there is a more complicated 
behavior at infinity, which however does not depend on n anymore. 

Before we state the RH problem that is satisfied by U, we first collect the 
properties of the functions gj that will be needed. 

6.2. The (7-functions. Observe that the functions gj are anti-derivatives 
of the Cauchy transforms Fj, i.e., 

(6.8) g'^{z)=F,{z). 

The asymptotic behavior of gj can be obtained by integrating the asymptotic 
behavior of Fj as given in Corollary 15.21 

Lemma 6.1. With the constant a as in Corollary \5.2\ we have that 

(6.9) gi{z) = log z + 0{z-^), 



(6.10) 52(2) 

(6.11) gs{z) 



'flogz- J^ + 0(z-i), Rez>0, 

Jlogz + ^-^^ 

'ilogz + g§ + 0(z-i), lmz>0, 

}logz + ^ + 0{z-^), lmz<0. 



|logz-^ + 0(z-i), Rez<0, lmz>0, 

|logz + ^-|^ + 0(z-i), Rez<0, lmz<0, 



as z ^> 00 uniformly. In (j6.9p ~ ()6.1ip every log z is defined with the principal 
branch. 

Proof. The asymptotics (j6.9p for gi follows from the fact that /xi is a sym- 
metric probability measure with compact support. 

The asymptotics (j6.10p - (|6.1ip for §2 and 33 can be found by integrating 
the asymptotic results for F2 and -F3 as given in Corollary [521 The constants 
of integration vanish, which can be seen from the fact that 



log(z - s)dfj,j{s) = log(2;)||/ij|| + / log(l - s/z)dnj{s) 

and log(l — s/z)dfij{s) -^ as z ^ 00. 

The additional constant ^ in (j6.10p in the third quadrant comes from 
the fact that we used the branch ()6.3p of the logarithm to define g2 which 
does not correspond to the principal branch of log 2; as z ^ 00 in the third 
quadrant. D 

6.3. The 0-functions. In the statement of the RH problem for U it turns 
out to be convenient to use certain functions (j)j. We recall that S,j is defined 
as the restriction of ^ to the sheet TZj of the Riemann surface. We define 

(6.12) </>j : C \ (M U ( - ioo, -ic] U [ic, icx))) -^ C, 
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as 



(6.13) Mz) = -l I (6(2/) -6(y)) dy, 



(6.14) M^) = I ' ' 



2/-!icfe(y)-e3(y)) dy, lm^<0, 



1 r^ 



(6.15) Mz) = j-ll fe(y) - e4(y)) dy. 

All paths of integration lie entirely (except for their starting values) in C \ 
(M U ( - ioo, -ic] U [ic, ioo)) . 

The importance of the (/)-functions is that jump properties of the function 
gj can be conveniently written in terms of the function (pj. We start with 
the jump properties of gi and the connection with (pi. 

Lemma 6.2. (a) We have for z G [—a, a], 

(6.16) ffi,+ (z) - 5i,„(z) = 20i,+ (z) = -2,/.i,_(z), 
and 

(6.17) 5i,+(^) +5i,-(^) - 52(^) - V{z) + ^r4/3|^|4/3 +^ 

^fo, ^e(0,a], 

|-27ri/3, ze[-a,0). 

(b) For z G M \ [—a, a], we /lei's 

(6.18) gi,+{z) - gi^^{z) = < 

I ZTTi, z E (— C!0, —a). 

(6.19) <7i,+ (^) +5i,-(^) - 52(^) - Viz) + ^r4/3|^|4/3 +^ 

_ \2(t)i{z), z G (a,oo), 

I 2(/)i^+(z) — Svri/S, zG(— oo,— a). 

Proof, (a) Let z G [—a, a]. Then we have 

(6.20) gi,+ iz) - gi-{z) = i / (arg+(2; - s) - arg_(z - s)) d/ii(s) 

= 27ri / d/ii(s). 
By ([OO]) . ([6l3]) . and the fact that 6,-(-z) = 6,+(-z) we obtain 

(6.21) <7i,+(^) - 51,- (^) = r (6,+(2/) - ei,-(y)) dy 



L 



z 
a 



(6,+(y) -6,+(2/)) dy 
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= 2</)i,+ (z) = -2(/)i,_(z) 

which proves (j6.16p . 

For z £ [—a, a] = S{fxi), it follows from the variational condition (|2.13p 
that 

(6.22) Re (^g,^^{z)+gi,^{z)-g2{z) - V{z) + ^r^/^\z\'/^ + A = 
and hence 

(6.23) 5i,+(^) + 9i,-{^) - 92{z) - V{z) + ^r4/3|^|4/3 + ^ 



(z-s) dfi2{s), 



= i / (arg^(z - s) + arg_(2 - s))d/ii(s) - i / arg 

J JiR 

where by the choice of branches in the definition of gi and g2, we have to 
use — vr < arg(2: — s) < vr in the first integral and — 7r/2 < arg(z — s) < 37r/2 
in the second integral. Clearly, 

(6.24) arg_,_(2; — s) + arg_ {z — s) = 

for z, s G R. From the symmetry of fi2 with respect to the real axis, we 
obtain 

(6.25) / arg(2; - s) d^2(s) = / (arg(z - s) + arg(z + s)) d/i2(s)- 

JiR Jo 

Since by (|6.3|) 

'0 if z > 0, 



(6.26) arg(2; - s) + arg(z + Sj - ^ ^^ if z < 
for s E (0,ioo), and since /i2(ilR+) = 1/3, we get from (j6.25p 

(6.27) 1^^ arg(z - s) d^2{s) = | ^J^g |[ ^ < JJ; 

Inserting ([OTll and (lOil) into (fOHD leads to dSHD. 

(b) Let z G M \ [—a, a]. Then we continue to have (j6.20p . and from this 
the relations (j6.18p immediately follow. By (j6.8p and (j5.3p we further have 

(6.28) 
^ (^g,^^(z) + 5i,_(z) - g2iz) - V{z) + ^r4/3|^|4/3 + ^^ = ^^(^.^ _ ^^(^). 

We obtain (I6.19P from (I6.13p . ()6.28p . and the fact that equality holds in 
()6.17p for z = iba. Note also that 0i_+(— a) = vri. D 

For the functions g2 and (j)2 we have in a similar way. 

Lemma 6.3. (a) For z £ S{a — /U2), we have 

(6.29) g2,+ {z) + g2,-iz) - gi{z) - g-iiz) = < -t ^ , ■ -i 

47ri/3, I] z £ (—100, — icj, 
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and for z G S{a — ^2) H iM±, 

(6.30) <7,,+(^)-52.-(.)±^r^/^kl^/^ 

= 202,+ (z) + 27ri/3 = -2(/.2,-(z) + 27ri/3. 
(b) For z G iM \ ^((T — ^2) = (— ic, ic), we have 

(6.31) <72,+ (^) -52,-(^) ± ^r4/3|z|4/3 = 2^i/3, z G iM±, 

and 
(6.32) 

/ N , / N ( ■. / N j'^Mz) ifz£ (0,ic), 

52,+ (^) + <?2,-(^)-5i(^)-53(^) = <„, . . . .,„ .. ^ / • n^ 

[202(-z) + 47ri/3 ^/zG(— ic, 0). 

Proof, (a) Let z G S{a — ^2)- Prom the variational condition (j2.15p it then 
follows that 

(6.33) Re {g2,+ iz) + g2,-{z) - gi{z) - 53(2)) = 0. 
Hence 

(6.34) g2A^) + 92,-{z) - gi{z) - g^{z) 

= i / (arg+(2; - s) + arg_(z - s)) dfi2{s) 

- i / (arg(z - s)) d(^i + ^3)(s). 

By symmetry we find 
(6.35) 

arg(2:-s)d(^i+/i3)(s) = / (arg(2; - s) + arg(z + s)) d(/ii+/i3)(s). 

Since arg(z — s) + arg(2; + s) = ibTr for z G iM-t and s G M+, and since 
(/Ui + ^3)(M+) = 27ri/3, we then obtain 

(6.36) / arg(z-s)d(^i + ^3)(s) = ±27r/3, z G iM±. 
Moreover, 

(6.37) / (arg+(z - s) + arg_ {z - s)) d/U2(s) = 27r/3, z G iM. 

Inserting ([06|) and (16371) into (lOij) leads to dOS]). 

By symmetry we have /U2(iIR+) = 1/3. Then by the definition of 52 and 
(j5.2ip we obtain 

/'°^ 97ri r^ 

d^2(y) = ~§~ ~ ^^W ^/^2(y) 
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(6.38) = ^ - / fe,+ (y) - 6,-(y)) dy - 27ri / da{y). 

"J JO JO 

By ([XT]) we have 

(6.39) 27ri /' d(7(y) = ±^^r^/3|z|^/^ for z G iM±. 

Jo 4 

For z G S'((T — /U2) n iM±, we also find since ^2 is analytic in (— ic, 0) U (0, ic) 

(e2,+(y)-e2,-(y))d2/= r (6,+(y)-6,-(y))dy 

J±ic 



(6,+(y) -C3,+(y))dy 

'±ic 

(6.40) =-2(/>2,+(z) = 202,-(z), 

where we used the definition (j6.14p of 02- Putting (j6.39p and (j6.40p into 
(16381) leads to (IO0|) . 

(b) Let 2; G (— ic, ic). Then as before 

(6.41) g2,+ iz)-g2,-iz)=2mj dfi2{y) = ^ - 2m J dfi2{y)- 
Now we have 

pz fZ 

/ d//2(y) = / dCT(y), 
JO Jo 

and so by (j6.39p 

(6.42) 2m /'d/i2(y) = ±^^r4/3|z|4/3^ fo^. ^ ^ (-ic,ic) n iM±. 

Jo 4 

Then (IHISTl) follows by (fOT]) and (IHIiSl) . 
Next, from (16.81) and (15.31) we obtain 



(6.43) _(52,+ (z)+52,-(^)-5i(2)-53(^)) 

= F2,+ (^) + i^2,- (^) - i^l (^) - F3 (2) = -6 (2) + CS (^) • 

Integrating (j6.43p and noting that the equality (j6.29p holds for z = itic, we 
obtain (I632D by (l6Tlj) . D 

Finally, for 5^3 and (p^ we have the following results. 

Lemma 6.4. We have for z G M, 

(6.44) 5 (^)+g3 (^)_5,(,) = /0' ^>0' 

I — 27ri/3, 2; < 0. 

and 

(6.45) 53,+(z)-53,_(.)^i^r^/3|,|i/3 

= 203,+ (^) = -203,- (^) + 27ri/3, /or z G M±. 
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Proof. From the variational condition (I2.17|) it follows that for z G M, 
(6.46) Reig3,+ iz) + 93,^iz) - g2iz)) = 0, 

and hence 



53,+ (^) + 53,- (z) - g2{z) =i I ( arg+(z - s) + arg_ {z - s)) d^sis) 

-i / arg(z-s) d/^2('S)- 



Combining this with ([Oi]) and (fOTl) gives (fOil) . 
Further, by symmetry we find /i3(M-|-) = 1/6 and so 

vri 

y ~ 

U,+ {z) gives 



2tt\ / d^3(s). 
Jo 



(6.47) g3,+{z) - 53,- (^) = 27ri / d^3(s) 

Now inserting (j5.22p and using ^3^_ (z) 

ff3,+ (^) - gsAz) = J-f (^3,+(y) - CsAy) - iVSr^^'lyl'^') dy 



TTl 

y 



'0 

\e3.4y)-e4,+(y))dy±^r^/^kr/=' 

4 



for z G M±. Combining this with ([613]) leads to ([05]) . D 

6.4. The RH problem for U. In the following lemma we state the RH 
problem that is satisfied by U. 

Lemma 6.5. We have that U satisfies the RH problem 

U is analytic in C \ (M U iM) , 

U+{z) = U-{z)Ju{z), zGMUiM, 

/I \ 



(6.48) 



[/(z) = (/ + 0(z-l/3)) 



z^/3 0/10 
10 \0 Aj 

Vo z-i/v 

as z ^ oo in the jth quadrant 

where the 3x3 matrices Aj are given in (|3.30p and ([3.3ip . The asymptotics 
for U are uniform as z —> oo in any region such that (|3.22p holds for some 
e>0. 

r/ie jump matrix Jjj is given by 

/g-2n<^l,+ ^ Q Q \ 

g-2n</.i,_ Q Q 

(6.49) Ju = „-2n<A3+ 1 ' o'^ 'S'(w), 



(6.50) 



Ju 



e-2"<^i> 


V 


e 




-2n</.3, 






+ 1 

e-2n03,-y 


/l e2'"^i.+(-) 

1 

e- 
\0 





-2n(/'3,+ 

e" 


\ 


1 

-2n.03,- J 


, ( 



onM\S'(/ii), 
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(6.51) 


Ju = 


/I 







g2n02,- 
1 





e2n</.2 


0\ 

+ 


, on S{a — fi2), 




Vo 





1/ 






/I 





0\ 


(6.52) 


Ju = 






1 

g-2n(/>2 




1 


oniM\5(o--//2), 




VO 





ij 


6.5. Proof of the 


jump matrices 


J(7 in 


the RH problem for U. 



Proof. From (j3.40p we see that the jump matrix for X on the real line can 
be decomposed into 2x2 blocks. The jump matrix (j6.7p for U can also be 
decomposed into 2x2 blocks, with two non-trivial diagonal blocks. The 
upper left block is given by 



(6.53) 



,n(sri,_(z)+£) 





X 





gn(92,-(z)-gi,-(z)) 



1 e 






.„(y(^)-|r4/3|^|V3) 

1 



e-'n-(92^+iz)~gi,+ (z)) 

V"(5i,+ W-Si,-W) gn(si,+(z)+si,_(2;)-92(z)-V'(z)+|r4/3|2|4/3+^y 







a"(9l,+(2)-9i,-(z)) 



which by the relations in Lemma 16.21 and the fact that n is a multiple of 
three, gives the equalities of the upper left blocks in (j6.49p and (j6.50p . 

For the lower right block in the jump matrix (j6.7p on the real line we find 
after a simple calculation 



(6.54) 
for z S 
(6.55) 



g-n(g3,+ (z)-g3,-(z)- 13^^4/3^4/3) 

e 

and 

-n(33,+ (z)-93,_(z)+i3^r4/3|2|4/3) 

e 



g"(93, + (2)+g3,-(2)-g2(2)) 
n(g3,+ (2)-g3,- (2)- i2^r4/324/3) 



g'^(93,+ {z)+g3,~{z)-g2 (z)) 

"(93,+ (2)-93,-(2)+i2^r4/3|2|4/3) 



for z E M_ . This proves the equality of the lower right blocks in (I6.49P and 
(j6.50p in view of Lemma 16. 4[ It is important again that n is a multiple of 
three. 

We finally come to the middle block in the jump matrix on the imaginary 
axis. By (j3.4ip and (|6.7p we find after a calculation 



(6.56) 



^(92,+(2)-92,-(2)+i^r4/3^4/3) 

nig2,+{z)+g2,-{z)-gi{z)-g3iz)) n{g2,+ iz)-g2,-{z)+'-^T^/''z*/^) 
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and 



for z E iM_. Then by Lemma 16.31 and the fact that n is a multiple of three, 
we indeed obtain (j6.5ip and (|6.52p . D 

6.6. Proof of the asymptotics for U. 

Proof. We deal with the asymptotic condition in the RH problem for U. 
Define matrices 



(6.58) 



(6.59) 



Hi 



H:> = — 




Ho 



Ha = — 




From Lemma |6. II and (j6.5p we obtain 

'1 

2^/3 



(^ 



\ 




/ + 



1 
H, 



,-2/3 



+ 0{z-') 



G{z) 



^+fj H,^' 



2/3 



+ Oiz 



-1- 



(6.60) Giz) = ^ Q ^,/3 Q 

V z^/^J 

for z — > cxD and in the jth quadrant with j = 1, 3, 4. In the second quadrant 

the situation is slightly different. Then we have 

(6.61) 

/z~i \ 

w^zVs 

LJzi/3 

V z^/^J 

for z ^ oo in the second quadrant. The extra factors a;^ in the (2, 2) entry 
and Lo in the (3, 3) entry are due to the extra term — f^ iii the asymptotic 
behavior of g2 in the second quadrant as given in (j6.10p . However these 
extra factors play no role in the asymptotic behavior of G{z)", since n is a 
multiple of three. 

Then by the asymptotics (j3.39p for X, 

(6.62) 

/I \ 

X{z)G{zr = {I + 0{z-^/'^)) 



1 

A, 



z-^/^ 

10 

\o zVsy 

1 



I + n 



H, 



2/3 + 0(z-l' 
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{I + 0{z-^/'')) 



/I 




Vo 



-1/3 






\ 




zV3 



/ + n 



A,H,A^' 



-2/3 



-(-')) GX 



as z 



(6.63) 



oo in the jth quadrant. A simple calculation then shows that 

01, J = 1,2,3,4, 



A^H^Af 




and hence 



(6.64) 



/I 



Vo 



-1/3 









1 





\ 




.l/3y 



{l + 0{z-'")) 



o> 



I + n 



1/3^ 



1 



/I 










AM^A:^ 



Vo 

as z ^ oo. Combining this with (|6.62p leads to 





Z-l/3 






-2/3 



\ 







,1/3 



(6.65) X{z)G{zY = (1 + 0(z-i/=^)) 



/I 





Vo 



-1/3 






\ 




,l/3y 



A, 



as z ^ oo in the the jth quadrant. Then the asymptotic condition in the 
RH problem (j6.48p for U follows by (j6.4p and the fact L is a diagonal matrix 
of the form (jG.Op . This completes the proof of Lemma 16.51 D 

Let us summarize the progress we made so far. Although the RH problem 
for U is not normalized at infinity, an important step is that the asymptotic 
behavior of U does not depend on n anymore. 

The jump structure for U is somewhat involved. Let us take a closer look 
at the 2x2 lower right blocks of the jump matrix Ju on M, that is, at the 
lower right blocks of ()6.49p and (I6.50p . We see that the diagonal entries 
are highly oscillating. Indeed, by ()6.15p . the fact that ^4^+(z) = ^3^_(2;) for 
z G S{a — ^2) and ()5.22ll we obtain 



(6.66) 



'/'3,+ (^) = -03,-(2) 



-vri/ dM3T '^r'lM'" 



for z E M, and hence (t)z,± is purely imaginary on M. 

An important step in the steepest analysis is the so-called opening of the 
lens. In this step jump matrices with highly oscillatory diagonal entries. 
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are turned into a constant (or otherwise nice behaved) jump matrix on the 
original contour and jump matrices with exponentially decaying off-diagonal 
entries on new contours on the lips of the lens. We do not go into details 
here, but let us mention that this construction works under the condition 

(6.67) ^Im(03,+ (x))<O 

which holds in our case by (|6.66p and the fact that ^3 is a measure with a 
density that is strictly positive on R. 

From (|6.5ip we see that the jump matrix for U on the the part S{a—^2) = 
(— ioo, — ic] U [ic, ioo) of the imaginary axis has highly oscillating terms on the 
diagonal. Indeed, by (|6.14p . the fact that S,z,+ {z) = S,2-{z) for z G S{a — ^2) 
and (|5.2ip we obtain 

(6.68) Cf>2,+ {Z) = -TTi f d(/i2-a) 

Jic 

for z £ S{a — ^2) H iM. In order to open the lenses successfully we now must 
have the condition 

(6.69) — Im</>2,+ (iy)>0, 

for y > c and y < c. The difference in sign in ()6.69p when compared to 
()6.67p is due to the different triangularity structure of the jump matrices. 
Indeed the middle 2x2 block in (|6.5ip is lower triangular, in contrast to the 
non-trivial blocks in (|6.49p and (j6.50p . which are upper triangular. From 
()6.68p we see that 

(6-70) —Im4>2,+ {iy) = -TT-—^{iy) + 7r——{iy) 

ay \az\ \az\ 

for y G M. And since iy G S{a — /i2) for y > c and y < ewe have that condi- 
tion (j6.69p is satisfied. Here we see the importance of the upper constraint 
II2 < (y. 

7. The third and fourth transformations U h^T h^ S 

7.1. Definition of tiie transformation U ^^ T. The next step in the 
Deift/Zhou steepest descent analysis is the opening of lenses. The aim of 
this step is to turn the oscillating diagonal entries in the jump matrices into 
exponentially small off-diagonal entries. 

We have to open a lens around each of the sets ^(^i), S{a—ii2) and ^(^3). 
The latter two are unbounded and we will treat the opening of lenses around 
these two sets in this section. The opening of the lens around the bounded 
set S{f^i) is more standard and it is deferred to the next section. 

The lens around S{a — 112) is opened as follows. In the discussion at 



the end of the last section we have seen that (16.69P holds for every y > c 
and y < —c. Then it follows from the Cauchy-Riemann equations that 
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Figure 7.1. The contours S2 and S3 enclose the inner re- 
gions ^2 and O3. Together with RUiM they form the contour 



Re (j)2iz) < for z in region around (— ioo, — ic) U (ic, ioo). In particular, one 
can show that, for some r > 0, it contains a cone 

(7.1) {z = x + iy gC I 0< |x| < 2r(|y| -c)} 

in its interior. 

Then we take the contour S2 around S{(t — ^2) as shown in Figure 17. 1^ 
so that 

• Re(/>2(2) < for z G 112 \ {±ic}, and 

• there exists an r > such that |x| > r{\y\ — c) for every z = x + iy £ 
S2 \ {±ic}. 

The lens around S{fi3) = M is opened as follows. Since ()6.67p holds for 
every x G M, we have that Re(/>3(z) > for z in a. region around M. The 
region is unbounded and it can be shown that it contains the cone 

{z = x + iy €C\0 <\y\ < 2r{\x\ + 1)} 

for some r > 0. We take the contour S3, as shown in Figure Ell such that 

• Re 03(2:) > for z G S3, and 

• there exists an r > so that \y\ > r(|x| + 1) for every z = x+iy G S3, 

see also Figure 17. li We may (and do) assume that S2 and S3 do not inter- 
sect. 

The contours S2 and S3 give rise to a partitioning of the complex plane 
as in Figure Ell The inner part of the lens around S{a — IJ.2) enclosed by 
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the contour S2 is denoted by Q2, and the inner part of the lens around 
enclosed by S3 is denoted by Q3. 

We define the 4x4 matrix valued function T by 



(7.2) Tiz) = U{z) 

(7.3) T{z) = U{z) 

(7.4) Tiz) = U{z) 

(7.5) T{z) = U{z) 

and 
(7.6) 



/I 0\ 

1 e2"'?^2(^) 

1 

\0 1/ 

/I 0\ 

1 _e2"</'2{^) 

1 

\0 1/ 

/I 0\ 

1 

1 

\0 -e-2"'?^3(^) 

/I 

1 

1 

\0 e-2"<^3(^) ly 



e ^.^..v^. ly 







z G O2, Rez > 0, 



z G ri2, Rez < 0, 



z € ils, Imz > 0, 



z G r^a, Imz < 0, 



r(z) = u{z) 



elsewhere. 



Then the 4x4 matrix valued function T is defined and analytic in C \ S^ 
where Tit is given by 



(7.7) 



St = M U iM U S2 U S 



3- 



7.2. RH problem for T. In the following lemma we state the RH problem 
that is satisfied by T. 

Lemma 7.1. T is the unique solution of the following RH problem 

T is analytic in C\ Sy, 

t+{z) = t_{z)Jt{z), zeJ:T, 

/I \ 



(7.8) 



r(z) = (/ + o(z-i/3)) 



1 

A, 



2^/3 

10 
\0 z-1/3 
OS 2; — > CO in the jth quadrant. 



V 



The matrices Aj are given by (I3.30p and (13.3ip . and i/ie asymptotic condi- 
tion in (|7.8p /loWs uniformly as z ^ 00 in each quadrant. 
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(7.9) 



The jump matrix Jt is given by 

/g-2n0i,+ 





1 



(7.10) 



(7.11) 



(7.12) 



(7.13) 



and 



(7.14) 



(7.15) 



(7.16) 



Jt 



J 



V 
/I 





1 

2n0i^, 






T 



Jt 



Jt 



Jt 





Vo 

/I 




vo 
/I 




/I 




Vo 







-1 



o\ 



1 

0/ 



o\ 



1 

0/ 



on 



on S{fi 



I), 



.\S{fii 











1 

I 0\ 
1 




o\ 




V 



on Sq 



V 



^2n(j)2 
1 










on, ^((T - /U2), 



on S3 



Jt 



Jt 



Jt 



/I 




Vo 

/I 




Vo 

/I 




Vo 






1 







1 

g-27102 




1 

g-2n(/>2 
g-2n(</,2+</.3) 






V 




1 





on (— ic, ic) \ ^3, 



o\ 




V 



on (0,ic) n Q3, 





1 

g-2n02 
_g-2n(02+(/'3) 



o\ 




1/ 



on 



-ic,0) n^a 



Proof. Each of the jump matrices (17.9l )- (j7.15p follows from straightforward 
calculations based on the definitions ()7.2p ~ ()7.5p and the jump matrices in 
the RH problem for U. Then the jump matrices (|7.9p - (|7.1ip are based on 
the factorization 



-2n4>3,- 
1 





g2n<^3,- 



1 

-2n4>3. 







1 

,-2n<t>3,+ 



UNIVERSALITY IN THE TWO MATRIX MODEL 53 

of the 2x2 lower right block in the jump matrix Jjj on M, see ()6.49p and 
(j6.50p . The jump matrices (j7.12p - (|7.13p are similarly based on the factor- 
ization 

1 e2'^'^2,+ j -\q I )\l Q )\{) 1 

of the 2x2 middle block in the jump matrix Jjj on S{(t — ^2), see (j6.5ip . 

The jump matrix (I7.14P is the same as the corresponding jump matrix 
(j6.52p in the RH problem for U . The jump matrix (jT.lSp follows from 
conjugating (16.521) with either (I7.4p or (I7.5p . 

The asymptotic condition in the RH problem for T follows from the def- 
inition (I7.2p - (l7.5p . the asymptotic condition in the RH problem for U in 
Lemma 16.51 and the fact that Re 02 < in ^2 and Re (^3 > in $^3 . For the 
latter facts, see also Lemmas 17.21 and 17.31 below. 

A more detailed analysis would show that the asymptotics for T actually 
holds uniformly up to the axes, in contrast to the asymptotics for U . To show 
this we would have to trace back the transformations Y^^X^^U^^T to 
see in particular what combination of Pearcey integrals is actually involved 
in the regions O2 and ^^3 near the axes. It turns out that the asymptotics of 
the relevant combinations of Pearcey integrals is uniformly valid up to the 
axes. We will not give details here. D 

7.3. Large n behavior of Jt- In the following lemmas we state result 
about the sign of Re 02 and Re 03 on various parts of the contour S-p. 

Lemma 7.2. For every neighborhoods U\c o-nd U-\c around \c and —\c there 
exists an 82 > such that 

(a) Re02(^) < -£21^1^/' for z G S2 \ C/±ie, 

(b) Re 02(2;) > £2 for z £ [-ic,ic] \ U±ic. 

Proof, (a) By definition of S2 we have Re02(-z) < for z €T,2\ {±ic}. Let 
us consider the behavior for z G S2 n I near infinity. Then 

(7.17) 02(z) = --j {Uy)-Uy)) dy=^ ^ z^/3 + 0(log(z)), 

as z — > 00 remaining in the first quadrant. Since Q2 contains the cone (|7.ip 
for some r > 0, we have that there exists a 6 > such that 

(7.18) 0<argz<iT/2-5 

for large enough z E S2 n /. Then we have from ()7.17p that 
(7.19) 

Re02(z) = -r4/3|z|4/3Re(^(u;-l)e4i'^'-g^/3) (l + 0{\z\-^/Hog\ 

as z — > 00 and z G S2 n /. By (j7.18p we have that 
(7.20) Re ({uj - l)e4iarg^/3\ = ^^^^ (-{avgz- ^72)") 
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< — Vonim sin — , 

- \^ 3 ' 2 y 

for z G S2 n / large enough. Hence there exists an £2 > such that 

(7.21) Re 02(2) < -£2^1^/^ 

for 2: G E2 n / large enough. By similar arguments in the other quadrants, 
we see that we can choose £2 such that (17.211) holds for z € S2 large enough. 
Finally, by continuity we can choose £2 small enough such that (j7.2ip holds 
for z £ T,2\ U±ic- This proves the first property 

(b) The second property follows from the variational conditions for ^2- 
From (j5.3p and (j2.16p we have 

(7.22) ReM^) =M'^92{z) - giiz) - gsiz)) 

= -iW^z) + t/^i(z) + [7^3 (z) > 

for z £ (— ic, ic). It is also clear that Re 02 is continuous on (— ic, ic). Hence 
the statement follows. D 

Lemma 7.3. There exists an £3 > such that Ke(j)3{z) > £3|zp^ for all 

Z£^3. 

Proof. The statement can be proved in the same way as we proved property 
(a) in the proof of Lemma |7.2[ D 

By Lemma 17.21 we see that the jump matrices Jt in (17.131) and (j7.14p 
converge pointwise to the identity matrix / at an exponential rate as n ^ 00. 
In fact, it shows that this convergence is uniform as long as we stay away 
from the points ibic. Since Re i;^3 > in Q3 we also have that (I7.15P and 
(|7.16p converge uniform to the identity matrix / at an exponential rate as 
n — ^ 00. 

By Lemma |7.3| we see that the jump matrix Jt in (j7.1ip converges uni- 
formly to the identity matrix I at an exponential rate as n — > 00. 

7.4. The fourth transformation T h^ S. In the next transformation we 
also open the lens Qi around S{fj.i) = [—a, a]. 

We use Si to denote the outer boundary of f^i so that Si consists of two 
contours from —a to a, one in the upper half-plane and one in the lower 
half-plane. Both contours are oriented from —a to a. See also Figure [7?2l 

Note that 

(7.23) (t>i,+ {z) = -vri / d^ui, z € [-a, a], 

J a 

Since we assume that /xi is regular we have 

d d 

(7.24) -— Im0i+(x)<O and -— Imc/fi _(a;) > 0, 
dx ' dx ' 
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for X G {—a, a), and we see that by the Cauchy-Riemann equations, there 
exists a region around {—a, a) so that Re0i(z) > for every z {—a, a) in 
that region. Hence the contours Si can be taken such that 

• Ke(j)i{z) > for z E Si \ {—a, a}, 

• El is contained in f^a. 

We define the 4x4 matrix valued function S by 
/ 1 0\ 

_g-2n0i(^) ^ Q Q 

10 

\ 1/ 

/I 0\ 



(7.25) S{z) = T{z) 



z e S^i, Imz > 0, 



(7.26) S{z) = T{z) 

and 

(7.27) 



10 

1/ 



z G ill, Iniz < 0, 



S{z) = T{z) elsewhere. 

Then S is defined and analytic in C \ S^ where 

S5 = M U iM U Si U S2 U S3, 

see Figure 17.21 for a sketch of S^. The next lemma gives the RH problem 
that is satisfied by 5. 

Lemma 7.4. S is the unique solution of the following RH problem 

S is analytic in C\ S^, 

S+{z) = S-{z)Js{z), 2; G S5, 

/I \ 



(7.28) 



5(z) = (I + 0(z-i/3)) 



^1/3 0/10 
10 [p Aj 

yo z'^/^J 

uniformly as z ^ 00 in the jth quadrant. 



The matrices Aj are given by (|3.30p and (j3.3ip and the jump matrices Js 
are given by 

/ 1 0\ 

-10 

1 
\ 0-10/ 

/I 0\ 



(7.29) 



Jsiz) 



z G S{pi), 



(7.30) 



Jsiz) 



e ^^M^) 1 
10 

1/ 



z G Si, 
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/ 


1 





0\ 



1 


(7.31) Js(z) = 






_g-2n(0i(2)+</.2(^)) 


1 

g-2n02(2) 


V- 


_g-2n((/>i(2)+</.2(2)+03{^)) 


g-2n(</.2(2)+<A3{^)) 


1/ 


for ze (0,ic)nJ7i, 








/ 


1 





0\ 


(7.32) Js{z) = 







1 

g-2n02(2) 




1 




K~ 


g-2n(<^i(2)+</-2(^)+<^3(^)) 


_g-2n{<?i2{2)+03{^)) 


l) 



for z £ (— ic, 0) n r^i, and 

(7.33) Js = Jt elsewhere. 

Proof. Since S differs from T only in a bounded domain, it is clear that the 
asymptotic behavior of S is the same as that of T. Hence the asymptotic 
condition in (j7.28p follows from the asymptotic condition in the RH problem 



The calculations that lead to the jump matrices are based on the factor- 
ization 

e-2"'^i.+ 1 \ _ / 1 0\ / 1\ / 1 

Q g-2n0i,_ j - |^g-2n0i,_ ^J \^_i qJ \^g-2n0i,+ i 

of the 2x2 upper left block of (j7.9p . We will not give further details here as 
this step in the RH steepest descent analysis is similar to the corresponding 
step in the RH analysis for orthogonal polynomials considered in [33]. D 



The following lemma deals with the sign of Re (pi at parts of the contour 

Lemma 7.5. Let U-a be a neighborhood of —a and Ua a neighborhood of a. 
Then there exists an ei > such that 

(a) ReMz) > ei for z G Si \ ([/_„ U Ua), 

(b) Re(t>i,+ {x) < -ei{\x\ + 1) forxeR\ {S{fii) U C/_„ U Ua)) . 

Proof, (a) This follows immediately from the continuity of (pi and the fact 
that Re^i > on S]^ -^ \ {0^,6^}. 

(b) Part (b) follows from the variational condition (|2.14p . Indeed, by 
(j6.19p we have that 

(7.34) Re0i,+ (x) = ^ Re (igiix) - g2{x) - V{x) + ^r^/^x^/^ + k^ 

for X G M \ 5(/ii). By ()2.14p and regularity of the measure we see that 
Re(pi{x) < for X G M \ S{fii). Moreover, for large x we have that the 
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Figure 7.2. The jump contour S5 in the RH problem for 
S in the one-cut case. 



dominant term at the right-hand side of (j7.34p is —V{x), which is a poly- 
nomial of even degree. Combining this with the continuity of Re 4>i we see 
that there exists an ei > such that 



(7.35) Re0i,+(x) < -ei{\x\ + 1) 

for all 2; gM\ 5(^1). 



D 



By Lemma 17.51 the matrix Jt in ()7.30p converges to the identity matrix 
I4 at an exponential rate as n — > 00. If one stays away from the endpoints 
of S{fii), the convergence is uniform. Combining Lemmas EH ED and [7]3] 
we also see that the jump matrix Jt given in ()7.3ip and (I7.32h converge 
uniformly to the identity matrix I4 at an exponential rate as n ^ 00. 

We also see that the (1,2) entry of Jt in (j7.10p converges to zero at an 
exponential rate as n — > 00. Again if we stay away from the endpoints of 
S{fj,i) the convergence is uniform. 



8. Construction of parametrices and the transformation S ^^ R 

8.1. The RH problem for M. If we ignore all exponentially small entries 
in the jump matrices Js in the RH problem for S, then we obtain the 
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following model RH problem for a 4 x 4 matrix valued function M. 



^.1) 



^ is analytic in C \ (M U S{(t - 
/ 1 ( 

4^) = M„(.) : 


, z E 


5(m), 


\ -1 oy 

/l 0^ 




\{z) = M^{z) 


10 
1 


; 


zGM\5(/ii), 


\o -1 o/ 
/l o\ 






V(z) = M_(z) 


0-10 
10 


5 


z G 5(cr - /i2) 


Vo [ 

\z) = {I + 0{z-^)) 


) 1/ 
/I 

zi/3 
1 


\ 






A 
lo A 






\0 





z-V3i 


\ -- 



uniformly as z — > oo in the jth quadrant. 

Here, the matrices Aj are given by (j3.30p and (j3.3ip . 

The solution of (j8.ip is not unique. To ensure uniqueness, we impose the 
additional conditions 

' M{z) = 0{{z T a)-^/^), z -^ ±a, 

M(z) =0((z^ic)-i/^), z^±ic. 



i.2) 



8.2. Construction of the outside parametrix M. The RH problem for 
M can be solved in the one-cut case by using a rational parametrization 
of the Riemann surface, which provides a conformal map to the Riemann 
sphere. 

We have S{fii) = [—a, a] and S{a — ^2) = iM \ (— ic, ic). The Riemann 
surface thus depends on two parameters a and c. It has genus zero. An 
explicit rational parametrization is given by the equation 



t^ 



sH^ 



(8.3) w + + 3 -., 

where the constants s > 0, t > are the unique positive solutions of the 
equations 



(8.4) 




-c. 



In this parametrization the branch points at z = ±a and z = ±ic correspond 
to w = ±s and w = ^it respectively. 
We introduce the following function 

(8.5) w ^ {{w'' - s^)iw^ + t^)f^\ 
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Figure 8.1. The cut r^r in the definition of the square root 
of the function w >-^ [{w'^ — s^)(w^ + t^)) 

defined and analytic in the complex tn-plane cut along 

(8.6) r^v := wi,+{Sifii)) U W2,- {Sia - /ia)) U «;3,+(5(/i3)). 

The square root is taken such that (|8.5|) behaves like w'^ as w ^ oo, and 
such that it changes sign when crossing the cuts. See also Figure 18. 1[ 



Proposition 8.1. The RH problem ()8.ip . ()8.2p for M has a unique solution. 

Proof. The uniqueness of the solution follows by standard arguments for 
uniqueness of RH problems, see e.g. |30J. We have to prove existence of a 
solution. 

Choose any basis Qk, k = 1, . . . ,4 of the space of polynomials of degree 
< 3, and define 



^.7) 



Nkiw) 



Qk{w) 



w;((w;2_s2)(w2 + ^2))i/2' 



Then each N^ is analytic in C \ Fjv and satisfies 

'" "' Nk,+iw) = -Nk-{w), weTN, 

^kiw) = 0(1), as t(; — > oo. 



9) 
10) 

11) 
12) 



Nk{w) = 0{w 



-1^ 



as w 



0, 



±s, 
'±it. 



Nk{w) = 0{{wTs)-^/^), asw- 
Nk{w) = 0{{w T ii)"^/^), as w 
Then the 4x4 matrix valued function z >-^ M{z) with entries 
.13) Mkjiz) = Nk{w,{z)), k,j = l,...,4, 
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where Wj, j = 1, ... ,4, are the mapping functions for ()8.3p . is analytic in 
C \ (M U S{a — /i2)), satisfies the jump conditions in the RH problem for M 
(due to the property (18. Sp of Nk), as well as the fourth root condition (j8.2p 
(due to (J8.1ip - (l8.12jl ). In addition we have as z 



oo, 



^.14) 



Mkj{z) 



f0(l), forj = l, 

\0(zV3), for J = 2, 3, 4, 

which follows from (l8:9]) - (l8lll and IKUh . 

Since the jump matrices have determinant one, it follows by standard 
arguments that z -^ det M{z) extends to an entire function. From ()8.14p it 
follows that detM(2;) = 0{z) as z ^ oo, so that 

(8.15) deiM{z)=az + P, 

for some constants a and /?. Since the polynomials Qk, 
linearly independent, a and P cannot both be zero. 
Define 



k 



1, ... ,4 are 



^.16) 



Aiz) 



/I 








\ 







^1/3 








(I 








1 





VO A 


\o 








,-l/3J 





for z in the jth quadrant. Then A is analytic in C \ (M U iM) and from the 
formulas (IS.SOp and (IS.Sip for the matrices Aj it easily follows that on MUiR 
it satisfies the jump conditions 



5.17) 



.18) 



A+{z) = A.{z) 



A+{z) = A.{z) 



(I 








0\ 





1 

















1 


Vo 





-1 


0/ 


/I 








0\ 








-1 








1 








vo 








1/ 



z G 



z G 



Comparing this with the jumps in the RH problem for M, we see that A 
and M satisfy the same jump conditions for z large enough. Then MA~^ is 



as z 



oo. 



analytic in a neighborhood of infinity. 

Due to JKT^ and ([8l6]) . we have M{z)A~'^{z) = 0{z^l'^) 
Therefore, by analyticity at infinity, we have 

(8.19) M{z)A-^{z)=C^O{z-^), z^oo, 

for a constant matrix C. Then det M{z) = det C+0{z^^), so that by (I8.15p . 
a = and /3 = detC. Since a and P are not both zero, we find detC ^ 0. 
Thus C~^ exists, and then it easily follows that 

-1' 



(8.20) 



M = C~'M 
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satisfies all conditions in the RH problem for M, as well as the fourth root 
condition I^M . 

This completes the construction of a solution of the RH problem for M 
M in the one-cut case. D 

8.3. Construction of local parametrices. The next step is the construc- 
tion of local parametrices near the branch points. Since we are in the one-cut 
regular case, the density of fii vanishes as a square root at the endpoints ita. 
As in the case of orthogonal polynomials [33^ ^^ the local parametrix will 
then be constructed with the help of Airy functions. Also for larger size RH 
problems Airy parametrices have been constructed, see e.g. [31 [HI [281 159] . 
The situation in the present case is similar, and so we will not give all details 
of the construction here. 

8.3.1. The model RH problem: Airy functions. Airy functions solve a model 
RH problem. Let yoi Hi and 2/2 be defined by 

(8.21) 2/o(s) = Ai(s), yi(s) = a;Ai(a;s), y2{s) = uj'^ Ai{uj'^ s) . 

where Ai is the Airy function and lo = e^'^^'^. Consider the 2x2 matrix 
valued function ^ 



(8.22) 


*(s) 


(8.23) 


^{s) 


(8.24) 


^{s) 


(8.25) 


^{s) 



yo{s) -y2{s) 
y'ois) -y'2{s) 
-yi{s) -y2{s 
-y'lis) -y'2{s 

-y2{s) yi{s) 
-2/2 (s) yi(s) 
yo{s) yi{s 
y'o(s) y'lis 



argsG (0,27r/3), 

args G (27r/3,7r), 
args € (— vr, — 27r/3), 
args G (-27r/3,0), 



Then ^ is analytic in the complex s-plane with a jump discontinuity along 
the rays arg s = 0, arg = ib27r/3 and arg s = vr. We equip these rays with an 
orientation as shown in Figure [8^ This figure also shows the jump matrices, 
which can be easily obtained from the definition (J8.22p - (j8.25p and the linear 
relation yo + 2/1 + 2/2 = 0. 

The asymptotic behavior of ^(s) is given by 

(8.26) 

as s — > cxD. 

8.3.2. The local parametrix near ±a. The local parametrix Pa is constructed 
in a small disk Ua around a. It should satisfy the jump conditions in the 
RH problem for S exactly and match with M on the boundary dUa- 
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Figure 8.2. The contour and the jump matrices in the RH 
problem for ^ 



Recall the function (pi from (j6.13p . which behaves near z = a as 



^.27) 



\Mz) = p{z - af^l + 0{z - a)) 



where p > is such that 
dm p 



(8.28) 
Then 

(8.29) 



i(x) = -{a-xy/\l + 0{{a-x)), x / a. 
ax vr 



/i(^) 



-l<Pi{z) 



2/3 



-ZG Ua, 



is a conformal map from Ua onto a neighborhood of the origin, so that 
fi[z) > for z G [/„ n [a,oo). We use the freedom we have in opening 
the lens around [—a, a] so that /i maps the part of Si n Ua in the upper 
half-plane into into the ray arg s = 2tt/2> and the part in the lower half-plane 
into the the ray args = —2tt/2>. 
Now define Pa by 

'^(n2/3/i(z)) o\ /e-"'^i(^)'^3 0^ 

/2y V h, 



(8.30) 



Pa{z) = Ea{z) 



for z G Ua, where a^ 

(8.31) 

for z G Ua- 



J*! ) and Ea is given by 



1-10 
Eaiz) = V^M{z) ( -i -i 
h, 



\0"3/4 



(n2/3/i(z))--- 



h 
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In the following lemma we prove that Pa exactly solves the jump condition 
in the RH problem for S in Ua and matches with the outside parametrix M 
on dUa- 

Lemma 8.2. The function Pa satisfies 

(8.32) 

Pa is analytic in Ua\ Jsi 

Pa+{z)=Pa_{z)Js{z), zeUaHJs, 

Pa{z) = M{z){l + 0(l/n)), as n ^ oo, uniformly for z £ dUa- 



Proof. The proof is standard. 



D 



The construction of the local parametrix P_a around —a can be done 
similarly. We can also use the symmetry to define P-a directly in terms of 
Pa as follows 



^.33) 



P-a{z) 



/I \ 




/I 


M 


0-100 
10 


Pai-Z) 


0-10 
1 






\0 -l) 




\^0 


-1/ 



8.3.3. The local parametrix near itic. The local parametrices near the branch 
cuts ibic are constructed in essentially the same way. We will focus on — ic. 
The jump matrices in the RH problem for S that are relevant near — ic 
(Trni) . (f7T3]) and (f7Ti]l . They are non-trivial only in the 2 x 2 middle 
block. 

We again construct the local parametrix by means of the Airy model RH 
problem and a conformal map that maps a disk U-ic around — ic onto a 
neighborhood of the origin. There is a small difference in the fact that the 
2x2 middle block in ()7.13p is upper triangular, whereas the jump matrix 
for ^ on the rays args = ib27r/3 is lower triangular. We can deal with the 
different triangularity structure by using 



^.34) 



$ = * 



1 

1 



instead of ^ in the construction of P-ic- 

The conformal map /2 is constructed out of (j)2, which behaves like 

3 



(8.35) 



-Mz) = P{Z + ic)=^/2g-3W4(i + 0(^ + ic))^ 



as z — > — ic. The fracional power is defined here with a branch cut along 
(— ioo,— ic] and so that 02 takes positive values on (— ic, 0). The number 
p > is such that (see also part (b) of Theorem (j2.3p ) 



^.36) 



d(q- - fi2) 
\dz\ 



(iy) 



-I 

IT 



y + c\'/\l + 0{y + c)), 



y 



-c. 
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Then /2 is defined on a small enough disk f/„ 

/o \2/3 

^•37) f,{z)=(-Mz) 



around — ic by 



with the 2/3-root taken so that /2(-z) > for z £ (— ic, 0) n U^ic- Then /2 
is a conformal map from U-ic onto a neighborhood of zero. We adjust the 
definition of the lens around (— ioo, — ic] so that the part of S2 n C/_ic in the 
left half-plane is mapped into the ray arg s = 27r/3 and the part in the right 
half-plane into the ray args = —2tt/3. 
Now we define P-ic by 

/IOO 

(8.38) P_ie(z) = E^-Uz) $(n2/3/2(z)) 

\0 1 

with <1> as in ()8.34p . and E'-ic given by 



1 











g-n02 (2)0-3 











1 



^.39) E^ 



^M[i 



/I 







0\ 
-i -i 
1-10 
1/ 



(n2/3/2(z))- 




The fractional power (/2(-z)) "'^' ^ is defined with a branch cut along (—ioo, 
and so that it is real and positive for z G (— ic, 0) n C/_ic. 

Then we have the following result, whose proof is again omitted. 



-ic 



Lemma 8.3. The function P_ 
(8.40) 



satisfies 



P_ 



is analytic in U-ic \ Js, 
,+ {z) = P_ic,_(z) Js(^), z G U-ic n Js, 



P-ic{z) = M(z)(l + 0(1/71)), as n —> 00 uniformly for z G dU^\c- 

The local parametrix P\c around ic can be constructed in a similar way. 
We can also use the symmetry in the problem and define P\c in terms of 
P_ 



^.41) 



as follows 










(I 





0\ 


) Pic{z) = 






-1 
1 








Vo 





-1/ 



p_ 



/I 





o\ 





-1 








1 





Vo 





-1/ 



8.4. The final transformation S ^^ R. Having constructed the outside 
parametrix M and the local parametrices P±a and P-ic, we define 

'm(z), zeC\{U±a^U±ic'J^s), 

P±a{z), ZGU±a\^S, 

^P±ic{z), zeU±ic\^s- 



iA2) 



P{z) 
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Figure 8.3. The jump contour S/j for R. Only the jumps 
on dU±a and dU±ic contribute to the asymptotic expansion 
for R. The jump matrices on the other contours are expo- 
nentially close to the identity matrix as n ^ oo. 



Recall that the 3x3 matrix-valued function (J3.14p - (j3.15p constructed out 
of the Pearcey integrals was also denoted by P. Since this function will not 
play a role anymore, we trust that the double use of the symbol P will not 
lead to any confusion. From now on P will always refer to the function 
defined in ^M)- 

Define the final transformation S ^^ Rhy 

(8.43) R = SP~^. 

Then R is defined and analytic in C \ (Y^s U dU±a U dU±ic), and has an 
analytic continuation to C \ Y^r, with 

(8.44) Yr = {dU±a U dU±ic U Si U Sa U S3 U (-ic, ic) U M) 

\([-a,a]Ut/±aUC/±ic), 

see also Figure 18.31 The circles dU±a and U±ic are oriented in the counter- 
clockwise direction. 

We obtain the following RH problem for R. 

Lemma 8.4. We have that R satisfies the RH problem 
R is analytic in C\ Tir, 



iA5) 



R+{z) = R-{z)Jr{z), ze^R, 

R{z) = I + 0{l/z), uniformly as z 



CXD, 
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where the jump matrix Jr is given on the various parts of Sjj as follows 

(8.46) Jr{z) = M{z)P±a{z)-\ z G at/±„ 

(8.47) Jr(z) = M(z)P±i,(z)-\ zG5C/±ie, 



(8.48) 

Jr(z) = M_(z) 



/l e2"<?^i(^) 0\ 

1 

1 

\0 -1 0/ 



M+(z)~\ zGM\([-a,a]u;7±a) 



^.49) Jfi(z) = M(z) 



.50) Jr(z) = M[z) 



^.51) Jr(z) = M(z) 



/I 0\ 

g-2n</.l{2) ^ Q 

10 

\ 1/ 

/I 0\ 

1 e2"<^2{^) 

1 

\0 1/ 

/I 0\ 

1 

1 

\0 e2"<^3{2) ij 

and on the various pieces of the segment (— ic, ic) we have 

(\ 0\ 

1 



M{z)-\ ze^l\U±a, 



M{z)-\ zGS2\?7±i, 



M{z)~\ ZGSg, 



(8.52) Jr{z) = M{z) 
for z £ (-ic, ic) \ {Qs U U±ic), 

(8.53) jRiz) = M{z) 

forzG (-ic,ic)n(Jl3\^i), 
(8.54) 



e-2"<?^2(^) 1 
\0 1/ 



M{z)-\ 



/I 0\ 

1 

e-2nM^) 1 

\ Q-2n{4>2iz)+4>3iz)) Q ]_/ 



M{zr\ 



Jr{z) = M{z) 
for z G (— ic, ic) n ili 



/I 0\ 

10 

_Q-2n{4ii{z)+4>2{z)) ^-2n4>2{z) j Q 

\ _e-2«{0l{z)+</'2{z)+<}!'3{z)) Q-2n(<f,2(z)+Mz)) Q l/ 



M(z)-i, 
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Proof. The formulas for J^ foUow by straightforward calculations. 

From the asymptotic condition in the RH problems for S and M in ()7.28p 
and (|8.1I) . respectively, we obtain 

R{z) = S(z)M{zy^ = 1 + 0(z-^/^) as 2 ^ cxD. 

Since the jump matrices Jr{z) are exponentially close to the identity matrix 
as z ^ oo (see (j8.56p . (I8.58J) . (I8.60p below) the improved error term R{z) = 
1 + 0(^-1) will follow. D 

From Lemmas 17.21 17.31 17.51 and the matching conditions for the local 
parametrices P±a and P±ic it follows that all jump matrices are close to the 
identity matrix if n is large. 

In fact, by (j8.46p and (j8.47p and the matching conditions in (|8.32p and 
(j8.40p it follows that there exists a constant Co > such that 

(8.55) ||j^(^)_/||<^, z Ci dU±ic^ dU±a. 

n 

Here we can use any matrix norm || • ||. The other jump matrices are ex- 
ponentially close to the identity matrix as n ^ oo. By Lemma 17.51 we have 
that 

(8.56) ||Jfi(z)-/|| <Ciexp(-nei(|z| + l)), z &^\{S{^,l)\JU±a) 
and 

(8.57) \\Jr{z) - I\\ < Ci exp(-nei), z G Si \ U±a 
for some constants Ci > and ei > 0. By Lemma 17.21 we have 

(8.58) ||jR(z)-I||<C72exp(-ne2(k| + l)), z G S2 \ C/±ic 

for some constants C2 > and £2 > 0. From the fact that Re(f)i{z) > for 
2; G r^i, Re(/!)3 > for 2; e r^a and by Lemma [72] we have 

(8.59) \\Jr{z) - I\\ < C2exp{-ne2), z £ (-ic,ic) \ U±ic 

where we can use the same constants C2 and £2 > 0. Finally, by Lemma [7. 31 
we have 

(8.60) WJniz) - I\\ < C3 exp(-n£3k|'/^), ^ G S3 

for some constants C3 > and £3 > 0. 
This leads to the following result. 

Proposition 8.5. There exists a constant C > such that for large enough 
n, 

(8.61) p(^)_/||<_^, z€C\Tr. 

Proof. The proposition follows from the above estimates on || J^(2;) — /|| from 
the arguments as used in [MJ Theorem 7.10]. D 

This concludes the steepest descent analysis of the RH problem for Y. 
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9. Proofs of Theorems 12.71 and 12.81 

In this section we prove our main results Theorems 12.71 and 12.81 We 
assume that /ui is one-cut regular. 

(n) (n) 

9.1. The kernel Kii • The kernel Kj^-^^ is expressed in terms of Y by for- 
mula (jl.30p . By following the sequence of transformations Y >-^ X i-^ U >-^ 

(n) 

T 1— > 5 H^ i? we obtain expressions for Kj^^ in terms of the solutions of the 
other RH problems. 



First transformation Y >-^ X . From the transformation Y >-^ X (|3.32p - 
(j3.34p . it follows that the right-hand side of (jl.30p can be written as 



(9.1) <^(x,y) 



27ri(x — y) 



(0 Wo^niy) Wi^niy) W2,n{y)) 



DnQ+in^/^Ty)e 



!;e,,(n3/v.))^-(^)"^-(^) 



ei,+ (n3/4rx)-ig+(n3/4rx)-iD^ 



-in-i 








for X > and y > 0. By arguments as in (|3.10p and (j3.19p we can rewrite 
(I9TD as 



(9.2) 



<^(x,2/) = ^-^^(0 e-"(^(«)-t-''^H''^) o)x-\y)X^{x) 



/1\ 






for X > and y > 0. Similar calculations show that (19. 2[) also holds for 
general x, y G M. 

Second transformation X h^ U. Applying (j6.4p and using (j6.5p and (|6.6p . 
we next rewrite (19.21) as 



(9.3) 



<\x,y) 



27ri(x — y) 



(o e""(^(^)-3^'^'l^l'^'-»i'+(2^)+»2{j/)-fci) o) 
X U^\y)U+{x) 



/A 






x,y e 
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Then by (16.170 and ()6.16p in case y G S{fii), and by ()6.19p in case y G 
M \ S{ni), we have 

Viy) - Ir^^V^' - giAv) + 92{y) -h^ -24>i,+ iy) + gi,+ {y) (mod27ri/3), 
Substituting this into (|9.3p . we arrive at 

:,n(gi,+ (x)-gi,+ (y)) 



(9.4) Ki'l\x,y) 



27ri(x — y) 



(0 e2'^'?^i.+(2^) O) 



xC/-i(y)[/+(x) 



/1\ 






, a;,y G 



Third transformation U ^^ T . The transformation (j7.4p acts only on the 
lower right 2x2 block. It does not affect the kernel (19. 4p . So we have 

(9.5) <)(x,y) 



, ^ p"(9i, + (^)-9i,+ (y)) , , 

We- -..^ - I (0 e2"'?^i.+(^) O) 



27ri(x — y) 



X T-\y)T+{x) 



/1\ 






, x,yGM. 



Fourth transformation T h^ S. Combining (|7.25p and (j7.27p with (j9.5p we 
obtain 



p"{9i,+ (^)-9i,+ (y)) 

(9.6) K^^I>{x,y) = ^^-^——^{-xi-a,a](y) e^"'^^'^^'^ « O) 



,(n). 



27ri(x — y) 



x5;i(y)5+(x) 



/ 1 \ 

X[-a,a](x)e2-<^l>+(-) 






V 



, x,yeR, 



where X[-a,a] is the characteristic function of the interval [—a, a]. 
Final transformation S ^^ R. Since S = RP, we obtain 

e"{9i,+ (^)-9i,+ (y)) 



(9.7) Ki'l\x,y) 



27ri(x — y) 



{-X[~a,a](y) e2"<^l.+ (?^) 0) 



X P^Hy)R^Hy)R^{x)P+{x) 



I 1 ^ 

X[-a,a](a;)e2"'^i-+(-) 


V / 



for x,y G M. This is the final formula. 
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9.2. Proofs of Theorems [221 and [2781 

Proof. Recall that R{z) = / + 0{l/n) as n — > cxd uniformly for z G C\ T,r, 
which we use in (|9.7p . Theorems 12.71 and 12.81 then follow by calculations that 
are similar to the ones in e.g. |19tl28j. 

Note that the factor e'^'^siM^yai.+iy)) in (pTH) disappears for x = y, and 
it also drops out of the determinantal formulas in ()2.20p and p.2ip . D 
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